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. INTRODUCTION 






In course of an attempt to apply direct vector methods to 
certain problems of Electricity and Hydrodynamics, it was 
felt that, at least as a matter of consistency, the foundations 
of Vector Analysis ought to be placed on а basis ind 
of any reference to cartesian coordinates and the main theorems 
of that Analysis established directly from first principles. 
‘The result of my work in this connection is embodied io the 
present paper and an attempt is made here to develop the 
 Difforvntial and Integral Calculus of Vectors from а point of 
view which is believed to be new. 
In order t» resline the special features of my presentation 
of the subject, it will be convenient to recall briefly the usual 
method of treatment. In any vector problem we are given 
certain relations among a number of vectors and we have 
dos deduce some other relations which these same vectors 
satisfy. Now what we do in the usual method is to resolve 
“each vector into three arbitrary components and thus rob it 
first entirely of its vectorial character. ‘The various characteris- 

vector operators like the gradient and carl are also subjected 
o the same process of dissection. We then work the whole 
lem out with our familiar scalar calculus, and when the 
‘analysis has been completed, we collect. our components 
| the result in vector language. [t is of course quite 
ae far ая it goes, the final vector expression of the result. 
iot only a succinct look ба? божо bat Ано а 
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suggestiveness of interpretation which they had been laeking |, 
in their bulky cartesian forme. Bat surely, strictly speaking, 

we should not call it Vector Analysis at all, but only Cartesian 

Analysis in vector language. In Vector Analysis proper we 7 
have, or ought to have, the vector physical maj 
our vectors represent, direct before our minds, and this charac- 
teristic advantage of being in direct close touch with the only 
relevant elements of our problem is sacrificed straight aw: 





tudes which. 














if we throw over our vectors at the very outset and work with 
cartesian components. We sacrifice in fact the very soul of 
Vector Analysis and what remains amounts praotically tom | 
system of abridged notation for cert - 
And operators ої cartesian caleulus which happen to recur avery 
now and then in physical applications. 
"The one great fact in favour of this plan ix that it affords і 
чи greater facility for working purposes, this facility по doubt 
arising solely from our previous exclusive fan ity with Carte- 
ian Analysis, But however nseful it might be in this direction, 











complicated formule 














and generally in making the existing body of Cartesian Analy. 
available for vector purposes, the process, I venture to think, 
is at best transitional, and the importance of the subject and 
the importance of our thinking of vector physical magnitudes 
direct as vectors, alike seem to demand that the whole of this 





branch of Analysis should be placed on an independent. basi 
But there iv а peculiar difficulty at the very outset. Histori- 
‘cally, most of the characteristic concepts of Vector Analysis, 
like the divergence and curl, had been arrived at by the physicist — 
and the mathematician in course of their work with the Cartesian t 
ealoulus and had even become quite familiar before the а 
~. possibility of Vector Analysis as a distinct branch of mathematics 
* by itself was explicitly recognised. ‘The vector analyst at 
й first then starts from these old concepts which happen also 5 
to be the most fundamental, bat it is his object right from the 
beginning to exhibit them no longer in their cartesian forms, 1 
but in terms of the characteristic physical or geometrical 
attributes which they stand for. Very often now a question 
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of selection arises from among the number of ways in whieh 

the same concept may be defined, different definitions being 

framed according to the different points of view from which 
% the subject is intended to be developed, The physicist—who, 
by the way, makes the greatest practical use of Vector Analysis 
aud- whose sole interest also in the subject is determined by 
the service it renders him in his work—aims, first of all, at 
his detinition representing most directly a familiar physical 
~ fact or idea; but, at the same time, and very naturally too, he 
holds the posibility of the definition yielding quite easily 
his useful working formule, of equally vital importance. But, 
unfortunately enough, these two distinct aims of the physicist 
are irreconcileable with each other, the most natural definition 
from the physieal point of view leads to the useful transfor- 
mation formulw of Physies only with the greatest difficulty, 
and the definition that yields these formulae with any facility 
is generally hopelessly artificial from the physical point of 
view." It is this irreconcileability of the two di 
of the physicist which, I venture t» suppose, 
sible for the persistence of cartesian calculus in Vector Analysis. 
For what is done is that definitions are first framed with « 
view to direct summing up of the simplest appropriate 
physical ideas, but then the necessity almost inevitably arises of 
seeking cartesian expressions for working purposes, for making 
Vector Analysis a serviceable and at the same time an easily 
manageable tool in the hands of the physicist 
- 1 may just illustrate my point by recalling how the usual 
nitions of the two most characteristic concepts of Vector 






































Mere to m paper by Mr. E. B. Wilson in the 
ER Amerikan Math Sos, vol. 16, on Unification of Vectorial 
where he criticises the artificiality in the detiuitions of divergence 
E Vy am lam wathematiciam, Bursl| Forti, which wero chosen 
‘solely with а view to thoir adaptability for establishing the working formula 

Jt Yomi Analyse wih ease. Thus Burali Forti's definition of divergence ie 
ea rp (a= ху), where a із any constant unit vector, 
rect connection with any intrinsic property of the 
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Analysis have been adopted from the simplest physical ideas 
which immediately identify them. Thus the idea of divergence | 
is taken directly from Hydrodynamies, and keeping before our 
minds the picture of fluid leaving (or entering) a small closed 4 
space, we define the divergence of a vector functior. at я point 

as the limit of the ratio, if one exists, of the surface-integral * 

of the function over a small closed space surrounding the point 

to the volume enclosed by the surface, a unique limit being 
supposed to be reached by the closed surface shrinking up to a 
point in any manner. 

Again, it ix found that some vector fields can be specified 
completely by the gradient of а scalar function, »o that the line 
integral fof the vector function along any closed curve 
(simply connected) space would vanish. "Thus the work done 
nil along any closed path in a conservative field of force, But 
in case the vector function cannot be жо specified, an expression 
of this negative quality of the function at a point ix naturally 
sought in its now non-evanescent line integra! along a small 
closed (plane) path surrounding the point. ‘The ratio of this 
line integral to the area enclosed by our path generally appton- 
ches а limit as the path rhrinks up to а point, independently of 
its original form and of the manner of it» shrinking, but dep- 
ending on the orientation of its plane. ‘The moreover has 
usually a maximum value, subject to the variation of this 
orientation, and а vector of magnitude equal to this maximum. 
value and drawn perpendicular to that aspect of the plane which 
gives us the maximum value is called the curl of the original 
vector function. 

Now these definitions, embodying, ах they do, the most — 
essential physical attributes of divergence and curl, must be - 
regarded as perbaps the most appropriate ones that could be 
given from the physicist’s point of view. But then comes the _ 





















































+ By the surface integral of a vector function, we always mean the м 
integral of it» normal component. 

+ By the line integral of а rector function, we always moan 
gral of ite tangential component. й 
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practical problem of deducing from these definitions the working 

rules of manipulation of these operators. The direct deduction 

being extremely difficult,* the already acquired facility 

% working cartesian calculus is naturally utilised for the purpose, 

and thus is reached the present position of Vector Analysis 
which T have already described. 

"The only way out of the dilemma would seem to be found 
by ignoring altogether both of these two specific interests of the 
physicist and looking straight, without any bias, to the req 
ments of Veetor Analysis as a branch of Pure Mathematics by 
itself. And paradoxical though й may sound, this course 
perhaps would ultimately best serve the phy: ends also. 
Atany mte, no free development of any science ін certainly 
posible, so long as we require itat every step to serve some 
narrow specifie end, 

We ask ourselves then, what should be the most natural 

i starting point of the Differential Calculus of Vectors? All our 
old familiar ideas of differential caleulus suggest at once that, 
whatever the ultimately fundamental concepts might be, we 
should begin by an examination of the relation between 
the differential of the vector funetion (of the position of а point 
P in space) corresponding to а «mall displacement of the 
point Р and thie displacement. This very straightforward line 
of enquiry I propose to conduet here, and it will be seen how 
їп а very natural sense we can look upon the diverge 
curl a» really the fundamental concepts of the Diferen 
“Caloulus of vectors, and how this new point of view materially 
Pais ples ойт analysis 
The first three sections are preliminary. In the first two T 
summarise the definitions of continuous functions and of Inte- 
штај» and briefly touch upon just those properties which I 
{ге in course of my work. The third is devoted to the 
Gradient | f a scalar function. ‘The real thesis of the paper I 

























































9 Compare, for instance, the аеону encountered by Mr. E. Cunning- 
im а peer ов the Theory Of Functions of a Heal Vector in the Proceedings 
ІА. oc, vol. 12, 1013. 















VECTOR CALCULUS 


begin in the fourth section where I consider the Linear Vector 
Function only with a view to developing what 1 have called the 
scalar and vector constants of the linear function, and although 
there is nothing very special about these ideas themselves, they 2 
will be found to lead very naturally to the concepts of Diver- 
gence and Curl and have been made here the foundati 
which my Differential Calenlus is built. The fifth section 
devoted to that Differential Caleulus and in the sixth 1 consider 
a few Integration Theorems and the divergence шиї - curl of an 
integral with a view to showing with what case these operations 
can be performed from my point of view. 

Notation, 

With regard to notation I use Gibbs’ here, although some 
of ite features are obviously meant to suggest easy ways of 
passing from Cartesian formule to vector, and eice verad, with 
which of course I am not at all concerned. 

For convenience of reference I reproduce the uotation for 
the multiplication of vectors. 

ПА, B are any two vector», 
the scalar product of A, Bix А.В= | A | | B [cos 6, and the 
vector product is А ж В which is a vector of magnitude | A I 
| В | sin б, and in direction perpendicular to both A and В; 
[A], | В| denoting the tensors of A and Band é the angle — 

| between them. 
Again, if A, B, C are апу three vectors, the notation — 
[ABC ] is used for any one of the three equal products 
A.BxC=B.C=xA=C.AxB=the volume of the parallelo- 
piped which has A, B, C, for conterminous edges- 
The following useful formula will occur very often : 
A x (B x C) = (A.C)B— (A.B)C. 
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I. 


CONTINUITY : DIFFERENTIATION OF A VECTOR 
FuNcTION OP л SCALAR VARIABLE. 


1. The functions we deal with will be mostly continuous, The 
position of w point P in «pace being specified as usual by the 
vector (= OP) drawn from a fixed origin О, the function /() 
{н said to be continuous at P, if corresponding to every arbitras 
rily chosen positive number 2, а positive number y (dependent 
on 8) ean be found such that 1/(r-+*)—/(r) «3, є bei 
vector satisfying the inequality 1<1 < y. The 
denotes the absolute value of the scalar if V is «scalar, and the 
tensor of У if V is a vector. 























ТЕ we construct the vector diagram as well, that is, if by 
taking another fixed point О" we daw the vector O'P' repre- 
senting the value of the vector function corresponding to every 
point the region in which the function iv defined, then Q 
being а point in the neighbourhood of P and Q’ the correspond- 
ing point in the vector diagram, our definition of continuity 
implies that any positive number й being first assigned, a 
positive number » ean be found «neh that so long ne the tensor 
of the vector PQ is less than у, the tensor of PQ’ will be lew 
than 8, Tt implies in other words that a sphere (of radius у) 
сап be described with centre P such that points Q’ in the veetor 
diagram corresponding to all pointe Q within (not o») f 
sphere will lie within a sphere of any arbitrarily small radius 
8 described with centre P”. 

We prove now that in the same ease the angle I"O'Q', 
that jx, the change in direction suffered by the vector fane- 
tion ean also be made arbitrarily «mall, For, in the triangle 
ora’, 
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Hence, xin P 





But PQ’ can be made arbitrarily «mall, and O'P’ is 
supposed to be finite. Hence sin P^O'Q* and therefore also the 
angle P^O'Q' can be made arbitrarily small. It follows that our 
direction as well. 

3. The funetion / (s) i» said to have а limit at P, if û 
being any point in the neighbourhood of P we have the same 
limiting value of the function no matter in what manner Q 
approaches P continuously 

If the fünetion is continuons at P, the limit exists at P and 
is equal to the value of the function at P, and conversely. 

Tf the limit does not exist at P, then either of two things 
may happen: (7) there may be different limiting values for 
different approaches to Р ; ог (i) there may be по definite 
limiting value for any approach or some approaches. In either 
сане the function is discontinuous at P. 

A third kind of discontinuity arises when the limit existe 
at P, but this limi! t equal to the value of the function at P. 

But, as has been remarked already, we shall concern our- 
selves practically with continuous functions alone, and an 
examination ої the sort of peculiarities we have just noticed, of 
what has been deseribed as the Pathology of Functions would 
‘be out of place here. The only discontinaity we shall соте 
across is the infinite discontinuity which arises when | / (r) | 
tends to infinity at P. 

8. Turning our attention then to continuous functione 
alone, we note that the sum and the scalar and vector products 






























‘of two continuous vector functions are continuous also. „The - 


case of sum is almost self evident, and we prove now that if Vy, 
V, are two continnous vector functions, the scalar 
pcia жй, . 




















CONTINUITY : DIFFERENTIATION. 9 


Let V’,, V’, denote the values of the functions ata point 
yc in the neighbourhood of the point м. We have only to show 
that for any positive number 8 assigned in advance, а positive 
number » can be found such that 


LVA VaV V, | <3, 
for all vectors є satisfying | «| <n. 
Now У РУ» -V,- У ЗУ E, - V HEY. 
HVV- V У, (У-Ү) 

+Уу.(У,'—У,)+(У,'—У,)(У„—У,), 
which is not greater than | У, | |У, -У, | + ГУ, | ТУУ, | 
ЖАРУ ЗУ | УУ, |. since the magnitude of the scalar 
product of two vectors is not greater than the product of their 
tensors. 

Hence, since the absolute magnitude of the sum of any 
number of quantities is not greater than the sum of their abso- 
Jute magnitudes, we have 

КУХ ЖУ: IIVI IV. 
з ГУ УУУ 

But since V,, V, are continuous, 

|У У, | «аву arbitrary 8,, provided only|«| < the 
corresponding m, and. | У-У. | < 

Of the two numbers у, wy, let a> 
thon provided | «| <m.. | У-У, | <8, and | У,'—У, | <3, 
and therefore | VV. VV, | < | Vi l8, + | Ve |8, +8,8,- 

Again, since | V, |, | У» | are supposed to be finite, given any 
positive number ô, we can always find 6, and ê, such that 5> 

ру, Та, + | Va lB +. 
EE; values now of 8, and 5,, we have 
лу, М | <3, whenever | «| <n, 
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4. Lf we consider in particular the continuous vector fune- 
tion of a scalar variable, we can easily adapt the argument of 
the usual scalar calculus and prove the theorems associated with 
continuity in that calculus. If r= (0) be the function consi. 
dered, ¢ being the scalar variable, we can prove in particular 
that if za (4) and ra=/ (/,) and p is any number 
such that | m, | <p< | э, |, then there is a value of ( lying 
between /, and /, for which | /(/) | =p. In other words, as ¢ 
Varie contincously from f, to /y, the tensor of м axsumos nt least 
‘once every value lying between ¢, and /,. 

Tt can further be proved that if F (+) is any continuous funo- 
tion, scalar or vector, of д" whero y itself i» a continuous function 
of а scalar variable £, then F is а continuous fanetion of 4. In 
case F is a scalar function, it follows that if F,, Fy 
the values of F respectively for (=f, aud /—/,, then as / 
varies continuously from /, to /,, F assumes at least once 
every value lying between F, and F, ; and when F is a vector 
function, it i» the tensor of F that assumes, ая f varies conti- 
muously from /, to fy, at least once every value lying between 
the tensors of F corresponding to "яз, and ¢=f,. 


5. Ш for the continuous vector function z—/(), а unique 
limit exists of /U)—/ as ¢ approaches 4 from either side, 



























(i.e, from values less than ¢ to £ and from values greater than 









® Prof, We have to show that if 8 i» assigned in ii ' 
found such that 
тової £8, 


mun ПРУТ де, є being the voias cl оваа bé a 
quaes 20) Мк кайнары канын З чання 
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4 with respect to £ and is denoted by : The function м in 
the same case is said to be differentiable at /. 

А function /) which is continuous and differentiable 
at all points in a certain region can in general be represented 
by'n curve in that region. The terminus of r will trace out 
the curve as / goes on varying continuously, and the vector 





М. will be at each point in the direction of the tangent to 


the curve аб that point, 





If F() is а continuous vector function of м, it follows -now 
from the last article, that the vector diagram of КО) eorres- 
ponding to points lying on any arbitrary but. continuous enrve 
P=) between any two specified points P and Q ix also a 
continuous curve lying between the corresponding points Р" 
and Q' in the vector diagram. 

б. Mem value theorem for rm). — ris a continuous 
and differentiable funetion of £ forall values of £ between any 
two specified numbers /, and /,, tl r, and r, being the 
values of r respoctively for /—/, and (=f,, we have 








mor = PUTA) 


where 0 is коте positive proper fraction. 


This iv proved, precisely as in the case of the corresponding 
‘theorem in scalar caleulus, by considering the function 
2 





=r 
2-7 (6—44) 


звис 

t=t, and (=(,, and of which therefore the differ- 
co-efficient will vanish at.some point between /, and ty 
oo is uunc 64), where 0 is а positive proper fraction. Т 





hich is continuous for all values of / between /, and /, and Ы 
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Graphically, if R denotes the vector to any point on the 
chord of the curve »=/(é) joining the points /, and /,, the 
‘equation of the chord is 


Rer, + FP (—/,); 


for obviously it represents for varying values / а straight line 
parallel to y, and gives Rr, at /—/, and Roy, at (= 
tys тоВа В”, say, represents then for any value 
of / the difference of the vectors to points on the curve and 
the chord corresponding to that value of ¢. If these vector 
differences are now drawn from the origin for all values of / 
from #, to f, their termini will give us another curve re- 
presented by R'—e(), which clearly is a continuous curve 
returning unto itself at the origin for (=f, and /=/,, and the 











vanishing of Ж for some intermediate value of / implies that 


in course of the journey of the terminus of R' from the origin 
and back to it again, there will be a position which will make 
the tensor of R' or »— R stationary. 





Qe 
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INTEGRALS. 


, T. The Vector Volume Integral—Given any finite con- 
“tinuous volume т, if for any convergent system of sub- 
divisions * ої the region, the vector sum = F,r,, where 
т. denotes any sub-region at any stage of the »ub-division 
and F, the value of F at any point within the sub-region 
ra, tends to a definite, unique limit as the sub-division 
advances, independent of the partieular convergent system of 
sub-divisions used and of the particular values of Е chosen 
within the sub-regions +., then this limit is called the volume 
integral of the vector function F through the volume r, and 















r 
je written f Fir. : 

м Without going into the question of the necessary minimum 

й condition for the integrability of Р, we prove without 





much trouble the only theorem we require in this eonnec- 
tion, vis, that if F is continuous at all points within a 
finite region, it integrable also through that region ро 
the continuity of F ensuring that if F., F.” are the values 
of Е at any two points in the sub-region r., the tensor of 
the difference of Е. and F,’ becomes arbitrarily small as 
the sub-division advances and each sub-region diminishes 
volume. 

A graphical representation of the vector volume integral 
may also be suggested here, Starting from any arbitrary 
point O', we lay down the vectors Е.т. as in the ordinary 
“polygon of vectors, In the t the polygon becomes a 
continuous curve, ending say in A.” Then the are O'A* will 

represent [| F|4* and the chord GA will represent our 
` volume integral Fir, Since F is supposed to be integrable, 







































Compare Hobson's Theory of Functions of a Real Variable, § 281 
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the chord O'A” will be unique, but we may bave an infinite 
‘number of curves like О"А" according to the different. orders in : 
which we may place the vectors Е.т. in forming the 

polygon. All these curves however will have the si 
length f] F | dr and the same terminal point A’. Further, to 








any point P' on any one of these curves there will corres- ۲ 
pond a unique point P in the volume т, and conversely, — * 
но that there ' 





n one-to-one correspondence between the points in 
the volume and the points on any particular eurve.* We 
have also 4= Fir, МОм is the vector O'P’, so that the 
tangent at any point P' on the curve is in the direetion | 
of P at the corresponding point in the volume. It may 
happen that f | F |r is infinite, but fF/r at the same time 
exists as a finite vector, Thus the curve may make an 
infinite number of conyolutions, but such that the terminal 
point А" is at a finite distance from O’. 1 
Given any continuous surface 
S in а region where the vector function F is defined, we 
form the scalar product Е. at each point of the surface, 
^» denoting the unit vector along the outward normal at 
any point to the surface, and the surface integral, in the 
usual sense, of the scalar function F.» over the surface we 
call the surface integral of vector function Е over the 
surface and denote it by f In other words, 
convergent system of sub-di 
is a sub-area at any stage of the sub-di і 
the value ої F.» at any point within the sub-area, tho unique 
limit to which zF,.».S, is assumed to fend as the sub- - 
division advances is called the surface integral of F- 

the surface S. But with the advance of the sub-d 
the arcas S, approximate to small plane areas 
tangent planes at points P, and tho vectors #8, u 












































в 5 
-volume and the pointe on а lime, for 
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may be regarded as representing these plane areas both in 
magnitude aud direction. We may replace therefore the 

% notation fFadS {by Еліт, do=ndS representing the ulti- 
mately plane element 78 both in direction and magnitude. 








Forming again the vector product of F and the veetor 
A,S, at each point and summing up for all points and 
passing to the limit in the same way, we have another 
surface integral fF xm/S or fF xls. This bus been called 
the skew or vector surface integral, ЈЕ. being the direct 
or scalar shall always mean Рае 
when we speak only of the surfs integral of F, referring 
to fF x do as the skew когіасе integral. 


If the surface S is supposed to be continuous and to pores 
morvover a continuous tangent plane at every point, the vector 
n would bé a continuous funetion over the surface, and if F 
1 issuppored to be continuons also, both F.n and F xn will be 
continuous functions and the scalar and. vector surface integrals 
of F over S will both exist. We shall always make thir 
supposition here. 

9. The tine integrals Given any continuous curve, if in 
any convergent system of subdivisions, p, i» the vector chord 
joining two consecutive points of division at any stage of sub- 
division in the system, and F, i» the value of F at any point 
` Pof the curve between these two points of division, then the 
‘unique limit to which SF..p, is assumed to tend as the sub- 

division advances is called the line integral of F along the 
curve AB. "The chord р. is obviously equal to the difference 
` n the values of 7 at the two points which it connects, and our 



















































д B 
integral may be denoted vy (Trae. It is further clear that 


` witli the advance of the subdivision, p, approaches in. direction 
‘the tangent to the curve at P, and if therefore we denote 
he unit vector along the tangent at any point of the curve by 
“the integral is the same as the line integral, in the usual 
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sense, of the scalar function F.t, and might be denoted by f F.t de, 
ds being the sealar element of arc. 


We might in the same way define the vector line integral 
J Fx dr, but this will rarely ocenr in the present paper. 

In any case we shall always suppose that the curve along 
which we integrate їх not only continuous, but also possesses 


n continuous tangent, so that t is a continuous vector function 
of the position of a point on the curve. 





The following properties of the line integral follow imme- 
diately from the definition. 





B Р в 
(a ( har F.dra | Par, P being any point on the 
A 4 Р 


curve AB. 

(ifi) Yt 1 is the length of the are AB and L, U tho lower 
and upper limits respectively of F.t for the curve AB (which 
limite are supposed to exist, though not necessarily to be 
attained), then 


B 
usf Е. dr SUI, 
4 











ying LSM SU, 





(i) Further, if M is some number sati 
в 

we mave f Fade MI ; and in ease F is continuous, so that F.t 
4 

in continuous also, the value M is attained by ЕА at some point 


в 
P of the curve, and we have {т = (Ру. 
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THE GRADIENT OF A SCALAR FUNCTION. 


10. Let F(r) be a continuous scalar function of the position 
of a point P (ОР=т) in a given region. If Q îs a point in the 
neighbourhood of P, such that PQ=ah where a is a nuit vector in 
direction PQ and h а small positive mber, then the value of 














Fat Qis F (гнав). If now the limit PL [F (r + ah) — 


F(r)] exists as a definite scalar function (different from zero) 
of n and r, mit would measure the rate of change in the 
value of the function for a displacement of P in the direction 
a. Supposing the limit to exist and denoting it by f (a, т), 
we have F(r-+ah)—F(r)=hf(a,r)+hy, where » and Б have 
the simultaneous limit zero. 












Now since h appears in the left hand side of this equation 
i only in the combination ah, and the first term on the right hand 
ide ix linear in h, it follows that this term i» linear in а also, 
‘The function f (а, r) then is a scalar function linear in a; it 
vanishes moreover with a, and therefore it must be of the form 
a. G (r), where G (r) is a vector function of г, independent 
ofa. 

If the limit in question exists now for every direction a. 
emerging from P, the rate of change of F(r) in any direotion. 
a is a.G(r), the maximum value of which obviously, for varyi 
directions a, is obtained when a is taken in the 
the magnitude ої the maximum value is equal to the tensor 
ої Ө. The vector G is called the gradient of the scalar function 
F. The gradient of а scalar function then may be generally 
defined as a vector in the direction of the most rapid rate of 
increase of the function and equal in magnitude to this most 
maid te 

D 
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11. Тһе same question may be looked at geometrically also. | 
We begin by proving* that if F (r) is а scalar fonction continu- 
Ous їп а certain region and docs not possess any maxima or 
minima jn the region, and if F is the value of the fonction at 





‘any point P of the region, then there passes through Pa surface 
оп every point of which F has the value РЬ і 


For, since P із neither а point of maximum nor minimum, 
1 alt the values of Е in the neighbourhood of Р cannot be greater 
than Fp , nor can all the values be less than F}, , and there would 


be points in the neighbourhood for which F is greater than Fp 


н" and there would be points also for which Е is less than Fp 


In the neighbourhood of P then, let Q be a point such that 


Fo>Fp sand В a point such that Fn <Fp « Now on account. 


t of the continuity of the function, a region can be constructed і 
^ about Q within which the fluctuation of the fonction is as small 
as we please. Hence there exist other points near Q for which 
also the value of the function is greater than Fp » Similarly 























d there exist points near R for which the function is less than Fy, 
i Hence the region consists of two distinct regions in every point 
| of one of which F>F,,, and in every point of the other 
К ЕР, 

l 


Again, since in passing from any Q to any R along а conti- 
ki nuous curve, F must on account of its continuity assume all 
intermediate values, it assumes the value Fp somewhere between. 


х: 
© This proof is adapted from the solution an example i Routh’s s 
VoL 11 (Bx. 2, $ 124), whore from the fact that gravitational 
“neither в maximum nor а minimum in freo space ін deduced 
| line cannot from part of а level surface. 
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є Qand R on that curve. Hence there is a continuous surface - 


of separation of the two regions at every point of which F=F 


V which proves our theorem. 





If now the surface possesses a tangent plane at P, we take 
A а point P" on the normal to the surface at P in its neighbourhood, 
Through this point P' also will pass a surface on every point 
on which F: / and PP’ will be normal to both the surfaces 





F= sb 





4 век, Supposü hat the limit. È 
pand F=F,,. Supposing now that the limit 


exists ая P’ moves continuously along the normal and approaches 
© P,a vector in the direction of this normal and equal in magui- 
tudo to the value of this limit is called the gradient of F() at Р. 
ГР" might be on the normal on either side of the surface, and 
it is assumed that the limit in question exist in either case and 

- tbat these two limits ave equal.) 

"To woo that the gradient so defined gives ux the most rapid 
rate of increase of the function both in magnitide and direction, 
we take a point Q in the neighbourhood of P on the surface 
on which P lio, Let Z P'PQ=0. Then the rate of increase 
‘of the function in direction PQ 




















establishes the identity of the definitions of gradient'in 
$ = article and the last- 


Wo denote the gradient by VF. If аР is the change in the 
i the position of P, we have 











— "where у and | êr | have the similtaneous limit zero. 

t If the shift êr is supposed to take place along a definite 
inuous c curve гом, then as we have seen ($4, p. 10) Е 
фе a continuous fonction of ¢ along that curve, and oar 


vet tel инен a be тана = f 
кыйы d ol - 
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- — Further, if /,, ¢, specify any two points К, L on the curve 
r= (0 and if Е,, Е, are the values of Fat К and L respeo- 
tively, we have by the Mean Value Theorem of $6, p.11. 


аР) 


ково (7, Mt 





че 


4 5 
where (Ж) yg tennetes the value of fF at some point М on the 


curve lying between K and L. That is to say, 


dr 


ds «м 





1—1) GP. ( 


In particular, if the curve is a straight line in the direction ої 
the (unit) vector 4 and 4 is the length of KL, so that KT: 
we may write F(r+oh)—F(r)=ha(VF).=ha. YF (r+a0h). я 
where © is a positive proper fraction ; or again, F (r-+a)—F(r) 
жалу (rh ба). 


We establish now the corresponding. integral formula, 
for which we prove first that if / (r) ім any vector function (not. 
necessarily continuous) integrable along a given curve AB, then 


4 р 











P being any variable point on that curve, tbe integral fro. че 
i a continuous funetion of the position of Р on the curve, 
Р 


neon the curve, we have. 
ч % 





—_— 
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a 
кока {л 
Р 


3 Q 
But (see p.16) | Sar | <UL, where U isthe upper limit 
x 


off for the curve AB and Z is the length of the aro PQ. 
Hence | F(r+«)—F(r) | &U/, 

and therefore | F(r4-*)— F(r) | can be made less than any 

ive number 4, Й only 2 i» жо chosen that &> UZ, 





arbitrary p 
or 1< Ê , which is always possible because U ін supposed to be 


finite. Again, since the curvo ix supposed to possess a continuous 
tangent at P (p.16), there is a finite portion of the curve about Р 1 
for which the arc measured from P and the corresponding: chord 
increase together. It is possible therefore to take a point Р" 











| on the curve in such a way that tbe are РР" «Чума also such 
ү that the ares corresponding to chords PQ which are less than 
| | FF | „are less than the are PP’ and less therefore than ү}. 
р It follows that for all vectors є satisfying ]«| « | PP^ |, 

where the arc PP'« Й, we have | F(r--) — F(7) | «5, which 


proves our theorem. GS £525 

14. We conelude that if /(r) is integrable along any curve 
in а continuous region, and if А is a fixed point and P any 
variable point in the region, then integrating along the various 
curves through A and P, we have any number of functions 


© ff dr, cach of wbich is continuous for a displacement of P on 
A 

















VECTOR CALCULUS 


the curve along which the integral is calculated in any case. 
Under certain conditions however (See $42), of which the ~ 
continuity of /' is one, the intergral is known to be independent 





n 
of the path of integration, and in this case therefore f. fudr 
к 


will define a unique continuous function F(r) of the position of 
P in space. Assuming these conditions to hold, and assuming 
in particular that /(r) is а continuous vector function, we shall 
prove here that /= VF. 

For, if / ix continuous, and since in accordance with the 
understanding in $ 9, the unit vector ¢ to any curve through 
A and Р is supposed to be continuous, /"/ is also a continuo 
function of the position of а point on this curve, Hen 
the point 7 and we consider another point 7^ on. the сигу 
neighbourhood of P, /" and ¢ being the values respectively of 
Sand (at r^, thon for аву arbitrary д, a positive number y can 
be found such that 

Ist -ft | <4, if p rr] <a = д 
whieh shows that | /^-/^ | lies between [| +8and [4| 
—68 ; every value, in other words, of /-¢ in the portion of the 
curve between r and У lies between [fot | +3 and | 74| —& _ 






‚ the integral J foor lies between 





If « denotes the vector r 


MC | E +8) and 161 — 8), 2 being the length of the 
are between г and 7. 


But F(r-9 — F() m fdr: 
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е) Е) =e So) | «| E, where 3 and | «| have the 
‘simultaneous limit zero. 
"This result, which is trae for all curves through A and Р and 
true therefore for vectors « drawn in all directions round РО 
shows that /= v F. 


If therefore we have any curve in the region, and лү, r, ane 
аву two points on the curve, we have 


Y» fs. 
f VEdro р уче Fry) — FO) 
n "n 
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15. The most general vector expression linear in ғ can 
contain terms only of three possible types, ғ, arb and esr, 
а, 6, є being constant unit vectors. Since r, avô and ex rare in 
general non coplanar, it follows from the theorem of the parallelo- 
piped of vectors that the most general linear vector expression 
can be written in the form 
Art ал+ ус xr 
where А, p, © are scalar constants, The constants p, у may 








moreover be incorporated into the constant vectors а and e and 
же write our general linear vector function in the form І 
Ф G)m Ar Raub ex г, P 


Y. only із а unit vector. 
Obviously, ф (г) is distributive ; 

that is, 6 (а+бу= Ф (a) +4 (В), = 

and further $ (kr) =k ¢ (r), where k is any constant. 

10. Theorem.—The sarface integral of ¢ (г) over any cloned 
P surface S bears a constant ratio to the volume T enclosed by | 

the surface, the constant depending only on the funetion but. 

being independent of the particular surface over which we 
integrate. 

"To prove this we ‘integrate separately the three terms of ф (г) 
over the epa 





в 
We kuow [Але =} дә =ЗАТ, 






Те енде fast, dr we break up the region 8 into 
cylinders with axes parallel to : 
cach of these cylinders like PQ will have an even number | 
diste with the surface, MEETS 













le 
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obvious as in that proof. If then de’ and de are the elements of 
surface on S enclosed by the cylinder PQ, we have 
6 b. do’=—b. de=area of the cross section of the cylinder PQ. 
Let P be the point r, then if = is the length of the 
cylinder, Q is the point r+z, 4 being a unit vector; and 
the sum of the contributions of de and da’ to the surface integral 
approximates, as the cross section of the cylinder diminishes, to 
„(++ ab) b. de’ +a. є b. do 
ie, to ab abdo? 
te, to abdr, where dT is the volume of the cylinder PQ. 
Hence tho whole surface integral 





s 
farb do=abT. 


Wo may just by the way note from the symmetey of the 
кена that farbdo=fbrado=wbT, and this result holds for 


Es any two arbitrary constant vectors a,b. That is for any two 


o Put 


constant arbitrary vectors a,b we have 

afrbdo=b far de=a'lT, which shows moreover that 
Jibda UT and fardo оТ. 

Generally therefore fradw=fardo=aT, a boing’ any 
constant vector. 

To return to our proof now, wo have to integrate fex rde. 
x В, ко that ¢xr=a.rB—B.ra (р. 0) 
“Hence fox rar ја. 800 — [В.о RT — a pt =0. 
We have therefore finally 
Solr) da — SNT Ha T 








ie, МЭ. а= Заа D, вау, 
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| Proof. We proved in the last article that foxrdo=0, o 
ig any constant vector. It follows that fer x de=0 (p. 6), 
'. fr X du 0, or frx do 0, because c is arbitrary, 
Also, Ја x do=b x fordo=bx aT [816] 
Again f(ex r)xdo=frede—ferde [p. 6). 
ScT—3eT=—2eT [§16]. 








s в 
Hence f(r) x де = Анас] x de 
—(axb+2)T=—CT, sny | 





с 





j $600) x d= * 300597). being a constant vector, (ИЙ 












our theorem i» established. 
| D and C which we find here associated with every linear 
veotor function, we shall always refer to as the scalar and 
vector constants respectively of the linear vector function. 
“18. We consider the fanction now 1 
Ar+bra—exr 
уко constant SA+ _ 





which їм immediately көеп to have the м 





constant ia—(axb+2e) which differs only in sign feam th 
wootor constant of (7). 

B are any two arbitrary месі 
[Mea Bae В) 

| — [Aat lan сха) m 8,4 (а), if we onll the new fonction 









"The two functions $(r) Ara rb ex г P 
and $()=\r+ rae xer 
тау on this account be called conjugate functions. 
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19, Since the scalar or vector constant of the sum (or 
the difference) of two linear vector functions is obviously the 
sum (or difference) of the sealar or vector constants of the 

stwo functions, it follows that the scalar constant of $(r)+ 
Фе) in 2D and its vector constant is zero; and that the 
seular constant of $(r)—4'() is zero and ita vector constant 
in 20. 

Obviously again the conjugate of (+) +¢'(r) in itself ; this 
funotion, that ія to say, ix self conjugate. And the conjugate 
of ФО) Ф) =F) ФС) = (ФС) Ф (е) ], which is the origi- 
тий function with the minns sign’ prefixed. Such a fanction 
has boon called skew or nnti-self-conjugato. 

‘Phe function фе) — Ф (9) in full is 

a.rb—b.ra+ 2e xr 
ax) xr+2exr (р. 6) 
=(a x B+ 20) x em C xr. 
Honce $(r) ean be written 
=O) 014 1090)—9(0] 
в O xr, where 26, (7) has been Written 
for the self conjugate function $(r) +4 СТ). 

Any linear vector function @(r) therefore can be expromod 
ns the sum of two other functions one of which is self conjagete, 
аж the same scalar constant as ФСГ) and no vector constant, 
d the other i» skew, has по scalar constant and the same 
| vector constant as (r). 
тһе resolt (7) 4 =C xr shows moreover that the vector 




















` ‘constant of all self conjugate functions is zero. 


20. "The vector constant of 4(7) may be exhibited in another 
manner, for which we ealeulate first the gradient of the scalar 
unction r.&(7). 
are) m Cbr) err) nC) 
=r +87) [$0) -$7)]—79 (0) 
наон) 487-98) 
mèn (+ C) en | Be | where y in 


i-i 
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а sealar number which has limit zero as | dr | tends to vanish, 
it follows that V[r.$(7)] —9(/)-- (ғ), and we can write — ^ 
$6)-1vIn$0)) 1C xr. 


Integrating now round any plane closed curve, we have © 9%” 


ХФ) r= Тов G))»dr MO x rdr. 
Since r.$() is single valued, the first integral on the right hand 
side is zero, because it is equal to the difference in the values 
of } гг) at the same point before and after circuiting. Гр. 
Also fr x Ar is twice the vector area enclosed by the cur 
a fact which becomes obvious by taking a new origin O’ in the 
plane of the curve. -For if OO =a, OP=r and O'P=p, we have 
r=pta and dr=dp, and frxdrafipt«)Xdp. Hence since 
[dp and therefore also fa x dp vanishes, the curve being closed, 
we have frxdr=fpxdp which is а vector normal to the plane 
of the curve and equal in magnitude to twice its area. 
Thus fo(r)dr=jfOxrdr=pCfrxdr=CwS, where S 
stands for the area enclosed by the curve and » a unit vector 
along the normal to the plane of the curve. 





























The ratio Н мне С.м does not then depend on the partie 


cular curve round which we integrate, but it depends on the 
orientation of the plane of the,eurve, on the vector я. This 
ratio obviously again attains ite maximum value when » is taken 
in the direction of C. The vector constant of @(r) then isa 
wector in the direction of the normal to that plane, round апу 
curve on which if we calculate the line integral of (7) the ratio 
of this integral to the area of the curve is maximum, and the 
magnitude of the vector constant is equal to this maximum 
ratio. 

21. There is just one bit of work more in connection with 
linear vector functions before we are ready for the Differential 
Caloulus of vector functions, | 

If a is any constant vector, a x (7) is of course also а linear 





vector function ої z. We proceed to find D, and C, the scalar. d 


and vector constants of a x 4(7). . 
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Integrating over any closed surface (eneloing volume T) we 
have by definition 4 
D, T — fa x фп) з 1 

fhlr) x der 4 
—a.O T, where C is the vector constant af Cr). " r 
Again, —0,T=f[a x $(7)) x do | 

$C )osle— fagl) (р. б), 

But fad(v)«lo—ofd(r).l0—aDT, D being the scalar constant 1 


of 4(2. 
Also галі is calculated immediately by breaking пр the 

volume into thin cylinder with axes parallel to o, as in $10, 
p.25. Thus if OP=r and PQ=xa, the sum of the contributions 
of the elements of area de and de’ at P and Q approximates, as 
the cross section of the cylinder PQ diminishes, to 

(reza) ado’ +¢(r)ade 
which again, since sado’ = — rade = vol. of the cylinder 

and Фла) =l) + 4(0), 
approximates to ф(а)дТ, T being the volume of the cylinder E. 
фіат 
Hence finally —C, T—$(«)T—aDT 
зо C m Da—$(a). 


| 
; 
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22. The Differential Caleulus of the scalar function of a 

| single (scalar) variable concerns itself with the rate of change of 
the fonction with respect to the variable. In considering in the 
same way the rate of change of a veetor function / (7) of the 
Position of a point in space, the first difficulty we meet with 
is that this rate of change is different for the different. directions 
in which the point P may be shifted. In fact, the position of 
Р being specified in the usual way by 
fixed origin, a change in the position of P of magnitude 4 
and in the direction of the unit vector a would be denoted by 
ай, and the change in the value of the function would be 
Jr ай) Ло). The rate of change then in the value of /at 
P for displacement of P in direction a is 


L f(r*ah)—f(r). 
4= Ly 


In Gibbs’ notation this is denoted by «V; we shall often 

denote it,—perhaps a little more expressively—also by df: 

In order that the limit may exist it is necessary that / should 

be continuous at І? in the direction а. For if / ix discontinuous * 
in this direction, then however small 4 might be, | /(r +a4)— 

Jr) | would be greater than a certain positive number 8 and 


therefore | 22-0 | can be made greater than any arbi- 





the vector r drawn from a 

















і trary positive number, and therefore the limit cannot exist. But 
| the continuity of f alone irection a cannot ensure the exis- 
ү tence of the limit, for which itis necessary that the fluctuation. 
Ж 


E 1 | енер | as a function of й should be arbitrarily 

small within a sufficiently small interval on the line а in the 

neighbourhood of Р. ‘The continuity therefore is a necessary — 

though not the sufficient condition for the existence of the edd 
question, М 
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t does exist as a definite 
n ofa and 7, itis clear as in $ 10, that that fuuetion 








will be linear in 4. We may denote therefore L }[/(r-+ ai) — 


k JF). by Ф (а, r), or more simply by ¢ (а), where (4) stands 
Tora linear vector function of a. The explicit presence of л in 
Фа, ғ) would serve to bring out the fact that the rates of 
change of /(r) are given by different linear vector functions at 
different points P in the region. 

If the limit exists for all directions о issuing from the point 
P,—for which it is of course necessary that /should be continuous 
at P—we write for any a 

daf or a. V f= ¢ (а), 
and (r+ ah) =/(r) + Ла) + hn, where 9 
| | has limit zero as А tends to vanish 

Or, since A 4 (a) $ (44) [ $15, p. 24), if f denotes the 
vector increment of / corresponding to the increment ди ofr, 

. буа) | êr |. 

23. If the shift бе be supposed to take place along a definite 
Continuous curve r=4(/), then we know ($ 4), that / would 
be a continuous function of £ along that curve and we would 
write from our relation of the last article that ~ 


P2] 

gant. 

“which, since $37) is a linear vector function of ди, 
М) C615 p.24) 


d 
і =). 

т, if £, ty specify any two points K and L on the curve 
f), and /, and /, are the values of / at К and Т, 'respec- 
then we have from the Mean value theorem of $6, 


fem 








a vector such that 
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where (A) denotes the value of ЧУ at some point M on the 
curve lying between K and L. In other words . 


fifi) [5]. 
Mf, of course, a defiuite ф exists at every point on the curve be- 
tween К and L. 

In particular, if the curve is а straight liue in the direction- 
of the (unit) vector a aud A is the length of KL so that KL= 
we may write /(r4-a4) Ло) ou, r4-a64) 
44./\е+ аб), 
assuming, of course, that d, / exists at all points on the line KL. 

24. We shall practically always confine ourselves to funo- 
tions which are not only continuous within a certain region, but. 
arealso such that 4,/ or (a) exists at every point P of 
the region for all directions « round that point. If we 
construct a sphere of unit radius with P as centre, then 
every point on this sphere will represent a detinite direction а 
issuing from P, and $ being known for P would mean that" 
corresponding to every point on the unit sphere we know the 
rate of change’ of / (at P) both in direction and magnitude. 
But as to the rate of change of / at P us a whole, we cannot as 
yet form any definite conception, not at least directly from our 
knowledge of the function $at P which only brings before our 
is а bewildering diversity of the rates of change for the 
` infinitely many directions round Р. What we naturally do there- 
fore is to have an idea of some sort of average value of ¢ (a) 
* for these directions a—average of $ (a) over the unit sphere 

round Р. 5 

"We consider then two kinds ої such an average value. 

‘Since а ів а unit vector, the magnitude of the component 
_ of ga) in direction a is a.#(0). We consider first the a 








































of this magnitude over the unit sphere, which is 
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where dS is the scalar element of area on the surface of the. 
sphere at the terminus of the vector а, and S is the whole 
„ surface 5. 


Since the radius of the sphere is unity, S is equal to 4= and 
їз therefore equal to ЗТ, where T is the volume of the sphere. 
Also adS=da the vector element of area, the normal to the 
sphere at the terminus of « being along а. Hence the average 
value 


= 1060-49 yp 


where D is the scalar constant of the linear vector function 
(a). [Dix а constant here in the sense of being independent 
of a, but is of course a function ог» and is different from point 
to point]. d 





Тһе average rate of change of / (r) then in the direotion of 
the displacement of P is proportional to the scalar constant of 
(a), and ti vernge rate therefore for each point of the sphere 
may be constructed geometrically by making the sphere bulge 
out uniformly outwanls from its centre P by an extra length 
proportional to D. 

‘To have an idea now of the average value of the tangential 
component of #(a) for all the points of the sphere, we consider 
naturally the average value over the sphere of the moment of 
Ф(а) about the centre P. This moment being ах ф(а), our 
average value 


леха ез) цс 


where C is the vector constant of (а), in the present instance 
of course the constant being a function ofr. The average 
moment therefore is in the direction of C and in magnitude is 
proportional to that of C. It follows that the average tangential 
component of ф(а) on the sphere,—the average of the component. 
that is to say, perpendicular to а,--ів perpendicular to C and 
` n magnitude is proportional to that of C. The vector constant 
of (а) affords us a knowledge of the average rate of change 
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of f(r) for any small displacement of r, perpendicular to that 
displacement. 

"The scalar and vector constants of $(4) therefore may be = 
regarded as supplying us with a basis of comparison of the rates 
of change of /(r) at the various points of the region 





in а sense the same purpose that is served by our old їп the 


ease of scalar function of a single (scalar) variable ;. 


35. The sealar D and the vector C being then so fundamen- 
tal in the Differential Calculus of Vector Functions, we hasten 
to exhibit them directly in terms of the function /(r) to which 
they belong, and we shall find incidentally how they are u 
mately identified with the well known Divergence and Curl of 
the vector function /(r). 

26. Divergence. We integrate. f /.de over any «mall closed 
surface surrounding the p P. At any point Q on this 
surface, PQ being ôr, the value of ^ 
еф) dm 


where A stands for | 8" | and » is a vector such that | 7 ' has а 
zero limit as À approaches zero. 


Hence f f. do — До (8r) 4m] de. 
But f fede =frfde=0, the surface being closed ; 


It). denm, D being the sealar constant of ¢ (êr) 
regarded as r vector funetion of де and T the volume 
enclosed by ааа 


Also Ја dee fh | т | dS, dS being the scalar magnitude of de 
< n JAds, where т is the greatest value of | 7 | - 


But fAd8—«T, where к is a finite number. 
Therefore, finde «дет. 














Hence Lffde—D «3^ 
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Now let the surface shrink up to the point Р in any manner. 
Then since when 4 approaches zero, all | 7 | 's and therefore 


also 7 tend to limit zero, it follows that L 2; fda exists and is 
equal to D. We have the following definition then— 

Enclose the point P by any small closed surface and calculate 
the integral f fdo. If the limit of HM Jda exists as the surface 


shrinks up to the point Р, which limit moreover is independent 
of the original surface and of the manner of ite approach to 
zero, then this limit is called the divergence of f(r) at P. This 
limit exists if g(a) exists at P, and in this case the divergence 
of /(r) is equal to the scalar constant of ¢(«), and may there- 
fore be taken (but for the constant factor }) as the measure of 
the average rate of change of f(r) corresponding to any dine 
placement of P in the direction of that displacement, the 
- average being taken for all directions round P. 

We shall always denote the divergence of /(r) by VA). 

37. It is perhaps possible for the divergence of a given 
function to exit at a given point P without $(2) necessarily 
existing there. Directly, the necessary and sufficient condition 
for the existence of the divergence at P is that within n suffi- 


8 
ciently small neighbourbood of P, the function n f fda should 


vary continuously for continuous variations of the surface S. 
In otber words, given any arbitrary positive num ber 5, й should. 
be possible to find a positive number у such that S,, 5, being 
any two closed surfaces round P and T,, T, the volumes 
enclosed by them, ] 
ГИК 
AS jde— pl jde 
ould be less than 8, whenever the surfaces S,, S, are entirely 


‘contained within the sphere with centre Р and radius о. It 
p 
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а matter for investigation now how far this condition necessarily 
Ц implies the existence of фа) at Р. We consider however just 
| now only those functions for which ф(а) exists at every point 








and for which therefore there is no question as to the existence > 
of the divergence. 





38. Curl.—Integrating ffxde over any closed surface in : 
the same way, we have 

if Sfx dam f [ fet (br) + hn) x do. 

+ But f fex do fex рй 0, the surface being closed i 


Jlr) x de — CT, C being the vector constant of ¢(3r) 
as a linear vector function of ёғ. 





Also fAqxde $ || ч| 48 E 
| «pas 
that is e Se T |. і 
D Sfxdo+C< пк, and therefore in the limit when tho 0 
surface shrinks up to the point P, we have 
1 * 
Lp Јуде С. 3l 


In general, if we find that the limit of Û wl чех exis as 





the surface shrinks up to the point P, the li 

boing independent of the original surface and of the sequence 
of forms taken by it during its approach to zero, then this limit 
is called the eur! of the function у) аё P. What we bave 
proved above shows therefore that if (a) exists at P, the curl 
does so too and in paie ma UU 
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29. The curl may also be exhibited in another (the more 
usual) manner corresponding to the property of the vector 
k- constant indicated in $ 20, p. 28. 

+ Thus, integrating f f. dr round any small closed plane curve 
3 surrounding P, we have 
Пе = Дурні бузку) di. 
But f fede fe- fdr=0, the curve being closed. 
e Solr) dr=C-nS, S being the area enclosed by the curve 
and n а unit vector normal to the curve. [$ 20]. 


Also, ЈА ае $ fh | n I ds, de being the magnitude of the vector 
jent of are dr ; 








_ aod |у | dr <7, JRAS, у being as before the greatest ||, 
and also fids=a8, where « is a finite number 


2 Фан у, 




















and we have 4 Sf-dr—O-n < 7a 
Hence, as the curve contracts to a point in any manner, 
I Pde approaches the unique limit Си. The limit moreover 


~ existe for all aspects of the plane area, for all vectors я. These 
various limits [for the different aspects of the plane have again а 
maximum value when я is in the direction of С, that is, when 
the plane is taken perpendicular to C and the magnitude of this 
“maximum value is equal to that of С. 

In case therefore Фа) exists for the function / (+) at a point 
(Р, it is indifferent whether we define the curl as we have already 
done it, or use the following definition — 

Having described any plane closed curve surrounding the 








int Р, if we find that the limit «1 fir, as the curve con- 





tracts to the point Р, exists and н, of the form of the 

ме and of the manner of its approach to zer», but dependent 
he ori of the plane of the curve, and if further, as 
n is varied, the various similar limits so obtained 
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for the different orientations all exist and acquire a maximum 
value for a certain orientation, then a vector drawn perpendi- 
cular to the particular plane which gives us the maximum value 
and equal in magnitnde to this maximum value is the curl of 
Jr) at Р. 

We shall denote the carl of /(7) by vx f(r) 

30. Some Transformation Formulae.—The explicit recogni- 
tion of the divergence and curl of a given function as the scalar 
and veetor constants respectively of the linear vector function 
9 (а, т), which defines the rate of change of the function for 
displacement in any direction 4, considerably facilitates the mani- 
pulation of these operators in practical work. This is what we 
proceed to illustrate. 

We shall in this article denote by w a continuous scalar 
function possessing a gradient at every point within the region 
considered, and U and F will stand for two continuous vector 
functions of which the rates of change at any point P will be 
же by the linear vector functions y (а, ғ) and ф (a,r) res- 
ely, so that the scalar and vector constants of ¥(a) will 
is у. and v x U respectively, and those of (а) will bey 7 
and V x У respectively, 

(9) То show now that 

VN) Y.V куму. 
aud V x(uV)euV XV + VuxV, 

Proof, Q being в point in the neighbourhood of P(PQ=ah), 
if the values of u and V at Q are w+é» and У+8У, the rate of 
change of «V for a displacement in direction а 


= 1 Ї (адн) (УВУ) -нУ) 
amo 

















1 
L 2 ру зкуванімву | 
RYE 3 


жна) ж Уа Ум. 2 
Bat V-(«V) and Welk are the scalar and vector 
regarded. 
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d therefore that the scalar constant of a.rb is a.b Гр. 25) and that ite 
vector constant is a x Б (р. 26), we have immediately 
E. VV) =u NR Ve 
) and V х(чУу=чу x V Tux V. 





(5i) To show that 
V(UxV)-V. V xU-U:g x V 
and V x(Ux V)e Uy. V—V V САУ) -ФИС) 
Proof. The rate of change of U x V for the displacement ah 
of the point P 


=L (0 +80) x (V BV) - U x v) 


=} [U xêV +80 x VY +êU xV] 


| =U x dla) +(e) x Vi 


апа we have to find the scalar and vector constants of this ax m 
linear vector function of a. 


Wo recall [$ 21, p. 29] that the scalar and vector constants of 
ax $(r) are—a:C and Ба—ф(а) respectively, Hence 


VaUxV)e—U: Vx VV. ухо 
and Y x(UxV)=U9 V—$(U)— V V U- (V). 
ug У-П УУ-УУ ОжЖУ vU. 
(ії) We may quite easily prove also the formula (ог V(U-)V 
. as given in Gibbs, Vector Analysis, p. 197, ris. 
VW V)=U x(V x V) Vx(V xU HU- VV-V-vU. 
Thus the rate of change of U-V for the displacement ah of the 
point Р 














=L} (+80). СУ НФУ) - U-V) 
=L ОзУ жи УВО ВУ) 


SUH + VG). 


— е 7 
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Now weknow that if this rate of change can be written in the б 


form «.G., then Ө = V(U- V.) 
But U фа)у+У Ма) Du. 

=U" [$(a)—$(2)]2- У [9(2) —V(0)] + Uda) + V (а) 
Р =U: (V хууха+У (V xU) x a-a'é(U) ea (V) Гр. 27) : 


=o (U x(V x V) Vx(V xU) eU) V)] 
Hence у (0 У) Ux (V x V)- V x (V х0) Ф007) МУ) 
Ux(VxV)-Vx(VxU)*U VV vU 





й We could write t 
: Since Офо) + Уа) could be written directly 

ca (U) x aw V), we have. 

VV) = GO) WC). 
In particular, V(a'V)-/(), a being a constant veetor. 


ame result in n more compact form. 














| A short note on Bilinear Vector functions. 
fle: 31. Before passing on to the Second Derivatives of the 
| vector function /(7), it would be necessary to consider v 
| briefly what are called the Bilinear Vector Functions. 
A vector function of two variable vectors linear in both 
called з bilinear vector function. iy 
1 Generally, a vector function of » variable vectors ene 
iù all of them is called an »-linear vector function. 





remains the same when the two vectors are interchan, 
rr denote the two variable vectors, the general s 








ra(r')A, where А is a constant vector and ¥ is 


WE RNC 
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= 32. Since the scalar constant of ab is «b (р. 25], the 

scalar constant of $7.7), regarded as а linear vector function 

? of ғ alone, ія БАЧ), or since the functions are all self- 
* conjugate, this scalar constant, 


-2EU40)—r).E90). 


of which again, regarded as a function of v, the gradient in 
m). 

It in now obvious а priori from the symmetry of $ (r,)—and. 
it is verified immediately also—that if we had calculated the 
scalar constant of фіто") regarded a» а function of " and then 
had obtained the gradient of this scalar constant with respect. 
to r, we would have got the same result (А). Hence, without 
ambiguity, we may refer to ЖА) as the gradient of the scalar 
constant of філо). We shall denote Sy(A) by Г. 




















38. Since again the vector constant of arb in axb (р. 96), 
he vector constant of (rn), regarded as a fune af v, in 
2907) ХА. We want to write down now the scalar and vector 
constants of this vector constant Z/(/)XA regarded as 
a function of ^. We recall that the scalar and vector constants 
of ахфіг) are—«C and Ра--ф(а) respectively, where D and C 
are the scalar and vector constants respectively of $(7). Hence, 
since the functions ФО) are all self conjugate and therefore 
their vector constants are zero, it follows that the scalar constant 
of БОМ) xA as a function of У is «ero. 


And its vector const 
" [Vs AL) +e (A) + eese] DD HAD, H. 
where D,...are the sealar constants respectively of ¥, (r), 
э. 
That is, the required vector constant 

* =7ф0)-ЖАР=г — SAD 
.— Amd here also the results would been the same if we 
` had enleulted the vector constant of $(,7), regarded as a func- 
оп of У, and then found the scalar and vector constants of 

vector constant as а function of 7. 



































© 
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We may say then without ambiguity that the scalar constant 
of the vector constant of ФСУ") Is zero and that the vector con: 
stant of the vector constant of $(r,^) = — ЖАП. We stall denote A 
this last by С". . 

34. For the bilinear vector function we have now to consider 5 
two conjugate functions. If we regard it as а function of r alone 
We have one conjugate function, and we have another when we 
regard it a» a function of ғ alone, We denote these two conju- 
gates by #,(r,7) and ¢',"(r,7’) respectively. These conjugate 
functions are not necessarily symmetrical Remembering that 
the conjugate of a 

И HW) = Ee GA Er 
КИСЛЕ EX) 
and 4G = ESO) 
and we propose now to seek for constants like [- and C' from 
these conjugate functions. 

Since a linear vector function (of one vector) and its 
conjugate have the same scalar constant and vector constants і 
| differing only in sign, it ix obvious that the sealar and vector 
| constants of $°, regarded a» a function of г and of ф/,' regarded 
| a» а function of 7 could be written down immediately from 
| the results we have already worked out for ф (7,77), but they 

would not obviously also furnish us with anything new. Also 
i Ф is only фо with r and r interchanged. We have to cal- — — 
| culate therefore only the scalar and vector constante of 4, 
regarded as a function of г. 
Since the functions V are all self conjugate, the vector 
constant in question is zero immediately, And the wealar 
constant = ZA: "осо SAD, the gradient ої which with respect | 
toris SAD. We denote SAD by г’. Thus г'= АР is gra- 
dient with respect to r of tbe scalar constant of 4', regarded 
(as «function of 7" and is also the gradient with respect to 
т of the scalar constant of ¢’, regarded as function Я 


















___ We thus have two independen: D the sym 
vector | 








ФО) meos 





© тези) 
EC "EAD 
à and С'=г—г 


35. Now we are in a position to consider repeated opera- 
tions of the derivative operators. 


Denoting L + Сеа) 707), the rate of change of /(r) 
4=0 


at P for a small displacement of P in direction а Uy œ, (а, г), 
. we consider first of all the rate of change of this function ф, (а, г) 
for any displacement of Р. If 0# is this new displacement, 
М af being а unit vector and 4’ a «mall positive number, then the > 
rate of change of Ф, (a, г) 





=L (sns rn) — s (ns п 










v=o 
Y Гү, Artak aA) оа) p Fran) еу 
=L (1 +“ эў: + 49 





#=0 4=0 hao 
=L + L FELCH mH + 0h) re у е ад) +A]: 


#=0 йно 
Assuming that a unique limit exists as 4, A’ approach zero, 
imi! vector function 
linear both in a and a’. We denote this bilinear vector function 
у фаба", а, r) and call it the second differential linear vector 
for f (r), $, (а, г) being the first. 
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By definition then 


dif =$, (а, г), 
du foa (a^, а, r). 








Tn the same way, 






4.4. 





y bpl (ew +08) Fe n) mf +a) +0). 
A-0420 
and this would be denoted by 4, (о, a’, 7). 
We see that М and d.d.’f differ only in the order in 
which A and 4' are made to approach zero and under certain 


circumstances, which may be investigated, tbe limit operations 
are commutative and then we should have 


dd f md d tf. 


case df, d.d. f. and did, 
n a finite region round P, the 
commutative property of the limit operations in question 
certainly holds and we have d, d. fed, d, f: 






We prove here only that 
all exist and are continuous wil 


For, applying the Mean Value Theorem of $23, P. 32 to 
the function 


f (r+ ¥ )—f (r) 
we have [CA +¥ ад) fr tahy] — [I(r tak) Ar) 
=1. Деко аву Лек абА)у), Ө being a positive 
proper fraction, 
Mul M [für a8 ¥ каві), applying the same Mean 
Value Theorem to /(7 +484), @' being some other positive proper 
fraction. 
In the same way we bave 
[Artak eal) Доза Ky] — Lr + ed) SO) 
каена M) еве" Py) 
Wd hd [rt a8 А+), 
8',, Ө, being also positive proper fraetious, 
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= Provided only then r+ah, r+ aJ! and r +a aM lie within 
the region where our conditions hold, we have it that whatever 
o. A, W might be, 
? Jr a eal) — Fr al) Artak) tAr) 
MILL f(r +¥ + авд) 
fre we" каб, Y). 

Since again d d f(r +40¥ адд) and du LIP e 
a0 ,A) are supposed to be continuous, they approach the same 
limit at r, and we have 4 fd, d, з or palaat) e da (ауа), 
Та other words, $, is a symmetrical bilinear vector function of 
ва" 

36. Since by definition $, (н'уу) regarded as a linear 
vector function of а" gives us the rates of change of ф, (ау) 
for the directions a’, the divergence and curl of this latter 
‘function аге the scalar and vector constants respectively of 
ф.о) as а linear function of a’, We proceed to show now 
how the second derived functions of /(r) can be obtained from 
d just as we obtained our first derived functions—the diver- 
gence and curl of /(r)—from $,. Since the divergence of /(r) 
а scalar, we сап have its gradient, and the curl being a vector, 
сап have if divergence and curl We consider these 
order. 

37. Let D stand for the divergence of f(r); D=. МО 


_ isa point in the neighbourhood of P, such that Рай, we have 
_ for the divergence of f їм Q. 

D+8D=V./+8/)= V.f Vf 
? 80 =.= (9, (ahr) +h) 


A 
a V DeL У Іа hr) 
«уреруя 6 































=; [5e] - X (or). фу being а 


function of є, 





r. e 
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| 

| 
f But У, (a,r) is the scalar constant of 9, ( 
linear function of a’ and is therefore equal to 
aY D=a.f , and since а is arbitrary, we have 
vpn-r, . 
where ГО is the gradient of the scalar constant of $,(a', 
38, In the same way, if C ix the curl of fat P and C+8C ін 





a,r) regarded ax a. 
Г 18 32]. Hence 





аг). 


the curl at Q, PQ as before being ah, we have 
€-v xf. 
0+80=V x (F8 ) 
овочу xpo x [4, (ahr) +h] 








Ll v x [$, (ahr) +A] 





хф, (а), 
which we know is the vector constant of фу Согог) regarded ак 
а linear vector function of a’. 
Hence V.C and V x C, which are the »calar and vootor oon- 
stants respectively of d, C as а linear function of a, are respectively 
the scalar and vector constants of the vector constant of $, (ar). 
The former we have proved to be zero [$ 33, р. 41]. 
that is, V. V xf is always zero; 
and the latter was found to be F — Г", that is 
Vx(Vxf)-vp-r' 
the gradient with respect to a of the scalar constant 
a.) regarded as a funtion of a’, We shall presently 

















б find an interpretation of Г ' directly in terms of f. P 
| 
39. Yt was seen [$ 30, (iii), p. 40] that if 
a df =$, (ar), 
then V(ej)—49,' (аг), where ф, (ау) is the linear vector » 
funotion of a conjugate to $, (a,r). 4 


Hence Y.Y (a: f у Умв, (ат) 
that is, Y.Y (a.f) in the scalar constant of 4.79, (ау) rogard- | 
ed as a function of a’. Ы сч» 







а 
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Now d,^ 





фара) фал) 
+ Therefore, B being an arbitrary constant vector, 


і Вал Qno) m Ll [Bh artak) Ahar) 
Mo 





J (ase (Bob) Br] 





рее Gee) e (1 














=а.ф,(а' Br). 
Hence 4,4, (4,r) in the conjugate of $, (o'a;r) with respect to 
а. That i 
dib, (ary m'a (a'r) 
Hence Y.Y (a.f "DS (834 р. 42) 





That 
tude of the component of f in any 
ponent of ГО" in that direction. 





the divergence of the gradient of the (scalar) magni- 
irvetion i» equal to the com- 








40. Generally, starting from any scalar function м, we may 
have two second derivatives, the devergence and curl of its gradi- 








ent G— Vw. 
l If G+8G devote the gradient at Q, РО ал, we have 
G48G= (u +) = TH+ GR) 
} £80 = vu) = Y [ha Vu +h] 
| =V[ha.G-+hn) 
; \ 2@,G=L1 9 (ha G-+hn)=9(a-G). 







‘The divergence and curl of G now are respectively the scalar 
and vector constants of 4.0 as а linear vector function of a. 


For shortness’ sake the divergence of О, that is, V.CVw) 
always devoted by "м. We prove now that the curl of 


ре 








к ОТИНДИ 
9 
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For denoting d,G by ф(о) for а moment, we know that — . 
V(a.G)=¢'(a) : and we have just shown that d,G— V(«.G) 





In other words, ф(а) is а self conjugate function and therefore 
its vector constant ix zero. 





Hence, V ж (Ум) =0. 


31. We have thus obtained two second derivatives of a scalar, 
funotion м, wis, У X Ум) and (Vu) or У и, of which the first 
vanishes for all functions w. a 











For the vector function f(r) we got three second deri 
СУ), V X(V xf) and YY Xf), of which the last vanishes 
for all functions f and the other. two correspond to two of the 
invariants [ and C' of the symmetrical bilinear vector function 


$a la’, а, r). 








But there was a third invariant Г of 4, (a^, a, r) which wo 

saw was related to f by the relation E 
У мат! or V*(af)m«[F '. 

«f ) ix maximum when a is taken in the direction 

ide of this maximum value i» equal to the 








showing that У" 
of p’ and the magni 
tensor of r’. This 
second derivative of the vector function /(r)—It i. 
which the direction is that along which if we caleulate the - 
б component of f(r), the divergence of the gradient of the (scalar) —— 
magnitude of this component is maximum, and of which the Д 
magbitude is this maximum value. This, as the relation 















уча) аг’ verifies immediately, i» of course that old. deri 4 
E vative which is so familiar to us in its Cartesian form — E 
5 8* , дб, 8") cerry. т 5 
1 (B+ 85735) (wit ej twk) $ 
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without consciously seeking for it, we arrivent it all the same 
from an unbiassed and straightforward examination of 


‚#9. 


"We obviously require now a new notation for this second 
derivative, for it is not deducible by any repetition of the first 
derivative operators, vie., the nt, divergence and curl, Ая 
all writers on Vector Analysis,—not excluding mathematicians 
like Silberstein (Vectorial Mechanics, Chap. 1) who constantly 
“advocate the exclusion of Cartesian decomposition from Vector 
Analysis,—have always defined it by its Cartesian expression 
which immediately suggest for it the notation v'*/(v), thi» 
eS) ix the notation that is invariably employed. With this 
“notation now we write 
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42. The characteristic properties of the divergence and curl 
lead almost immediately to the integration theorems 


в T 
дф ffde—[pv4de, 


в T 
(DD [fxde—— [уху , 

i g over Ше surface S and through the 
closed region ; and again 
=f хуло. 
the surface integral here extending over the surface of any 
finite unelosed region and the line integral round the contour of. 
the unclosed surface. ‘The function /(r) in all cases is supposed 
to be finite, single-valued and continuous at all points of the 
region of integration. 

‘To prove the first theorem :—For any sub-division of the 

region into smaller closed volumes, we know by the usual 
argument of the integrals cancelling over the interfaces, (the 








- continuity of / ensuring the equality of the values of / at cor- 


responding points on the two sides of an interface) that 





8. ran В, 
[f4e— A f fae 
8, denoting the Se surface of any one of the sub-regi 


less than any arbitrary number 4, then 
fae p. cox [ $26.) 
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8 
J fde=SD, +, + 50.7, 


+ „ _ Det us pass to the limit now as the sub-division advances and 
k^ й therefore diminishes indefinitely. 


p Now Бут, is always less than з Sr,, 9 being the greatest 9, 
5 вт. «9T. 
But T 





upposed to be finite and in the limit у =0. 


7e dn the limit Sy, r, =0. 





Also, by definition, the limit of SD, r, inf Ddr or To fac. 


в H 
Hence f fdo= ГУ far. 


In precisely the same way we prove that. 


8 А 

J frdo=— [у хуй. 

То prove the third theorem, we break up, as usual, the un- 
closed surface by а network of closed curves. ‘Then having а 
dofinite convention as to the sense in which the line integrals are 
to be calculated round these curves, we prove first of all in the 
usual way that the line integral round the original contour is 
equal to the sum of the line integrals round the closed curves 
that have been drawn on the surface. Hence, using the same 
sort of argument as used above for proving the first theorem, 
and by a reference to §29, we prove quite easily that 

Sfar=f V fo. 
Corollary 1. For a closed surface fV xf4e—0, and this 
another proof of the theorem Y.Y xf=0 [§ 38, p. 48); 

r the closed surface may be taken as small as we please, and 
_ then integration theorem () will prove the result. 

"The line integral round every closed curve in 















у вний 
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of the region. But if the line integral round any closed curve 
like ADPEA drawn through the two points A and P is zero, 


Р 
that means that the line. ا‎ flr are the same, whether * 


we use the path ADP or the path AEP. Similarly 
connect A and P by any other path like AFP, then since 
the line integral round the closed curve ADPFA is zero also, 
it follows that the line integral for the path AFP is equal 
tothat for the path ADP. We conclude that the condition 





Р 
(promised in $14, p. 22) that fdr may be independent of the 


path of integration and may therefore define a unique. function 
of P, (A being а fixed point) is that the curl of / should be zero 
at every point of the region considered. 

43. 1f we put /=ou in the theorem (i) of the last article, 
а being any constant vector and м a continuous scalar funetion 
possessing a gradient at every point of the region considered, 
we have 








[nado f У «анудо 
a,u de ; (8 30, (0) p. 38] 
that is, for any arbitrary constant vector s, we М 
fido a. fS udr. 
| Tt follows that fude= fV udr. 
| Again putting /=a in the theorem (вії) of the last article, 
we have 








Jua. dr [V x (wa), de 
fv«xede — [$30 091, 
Ja. dex V», 

the integrals extending round the contour and over tbe surface 
respectively cf any unclosed surface. 

Since а is a constant vector, we may write уз 1 
a. fudr=a. fdox Yu 
ранах Ju, a being arbitrary. 
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44. Differential of an integral—Sivce the calculation of 
the gradient, divergence and curl all depend on that of the 
differential of the function considered, it is necessary to have 
formule for the differentials of functions given in the form of 
integrals, before we can get the result of operation on them by 
these derivative operators. 


Let / (r, А) be a vector or scalar function, A being an 
arbitrary (vector) parameter of the function. 


Consider the volume integral f,/ (r, А) dr ої the function, 
the region of integration being bounded by the surface 
F(r, А) 9 0, where F is а scalar function. The integral of course 


is a function of А alone, say x(A)- 


Imagine now а «mall increment 8А to be given to A and 
lot rand М denote the volumes bounded by F(r,A)=0 and 
F(r, А 8) 0 respectively. Then 


о) f (rede fl f (r, л) 





ле ABRIL f (0s ۸+8۸) de—f f (rA) de 


Е [f (A+B) ©. A) нні, f Cr Ne Bs 





Di Jr, Х-+8Х)4+ denoting that this integral is to be taken in 





the region between the surfaces F(r, А)=0 and F(r, А4 8А) =0, 
d + being reckoned positive or negative according as it is outside 
or inside the surface F(7, А) 0. 

Now as | 8х | becomes smaller, the surface F(r, 4+6\)=0 
a F(s,4)=0 and the volume between them tends to 
become a thin shell distributed over this last surface; and 
if rand r+8r аге corresponding points on the surfaces 
Fr) =0 and F(rA +5۸) =0, the volume dr of this sbell resting 
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On the element of area de at ғ on the surface F(r,A)=0 
tends to 


rade or r.ndS * 


Hiat in, to TE as, 





m the anit vector along the norma! to F(r,X) 0 being Y-F. 
z TV.FT' 


if Y.F denotes the gradient of F(r, А) regarded as a function 
of r alone. 


Hence, as | 8А | diminishes, 


f уез A+A) de 


в змов 
n tends to f f(r, A+B) SONIS dS, the surface integral being 


taken over the surface S of F(r, X) =0 
То express this integral now explicitly in terms of ёл, we note 
generally that 

F(r+êr, A+8A)— ВА) 

=F (r+ 8r, АВА) Ее, A+) + F(7, A+ 8А) ЕСА) 

me, F(r +r, АВА) ВХ АКТУ AHOA) [$ 19, p. 20.) 
which, by sufficiently diminishing |A | and | ër] may be 
made to approximate, to any arbitrary degree of accuracy, to 

êr. V. Fir, А) + ВА УА Flr. A). 
if of course both VF and VAF are supposed to be continuous 
functions ої r and A, . 

If now r+8r be supposed to Ье a point on the surface 
X434, in the neighbourhood of the point г on the surface A, 
then 
| F(r+êr, А4 8А) =0 and F(r, А)=0, 
and therefore Sc. „Е (rA) + 5А. АЕ (У.А) oan be made arbitrarily 
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Hence the integral [d f Cr Ae BX)de further approximates to 


че A+) 72 


vr 


the degree ої approximation remaining the same. 





We write finally therefore h 


êx(\)=8 f f (r, X) de 
= f (f(r, АВА) ў (r, A)) de 


$ AF. 8A 
= f e.24 Ts rp 
our assumptions about the function F being that both т „Е and 
YaF are continuous functions in each ої the vectors r and А, 
45. Suppose now / to be a continuous scalar function 
possessing a gradient at every point within the region of into- 
gration. "Then by the Mean Value Theorem of $ 12, р. 20, 













a ен [A VAGA 0,83) Me 


осела дунд, Dy PAE ав. 





6, being a positive proper fraction. 
“If we further assume the continuity of Va / in А, then since 
. the surface S and the volume T are both supposed to be finite, 
the difference between the right hand side of the last equation and 





8 
Ут fva fe хо fir. X) var ozs, 


Rl vui f SE, 1481 
xero as | 3 | a Hence 


Е. - М 
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A6. Suppose next that / is a continuous vector function. 
Let def denote the rate of change of /, regarded as a function 
of A, for an inerement of A in the direction of the unit 
vector a. y 

If we assume now that da fis a continuous function of A, 
then using the Mean: Value Theorem of $ 23, p. 32, and 
remembering that the surface S and volumeT are finite, we 
have, by an argument similar to that of the last article, 





dax( A) = fdafile— f f NS 


But the divergence and curl of ХОХ) are the scalar and 
vector constants respectively of dax(A) regarded as a linear 
vector function of a. 

‘Now since the scalar constant of the sum of any number of 
linear fonctions is the sum of the scalar constants of those 
functione, and since further the scalar constant of de / іє 
V f, we have the scalar constant of fda /4 T —fVA'// T. Henoe, 
[remembering that the scalar constant of arb=a'4 (p. 25)). 





"c d 
vxo) vif fi хмтафую і / B et 


© Similarly, 


E Ж. в 
дк Повне VAX іно چ‎ t 
| 47. IfA is not involved in the equation of the bounding 
surface, and its variation therefore does not affect the region of 
integration, we bave simply, for a scalar function /(r) 


y vy ffir = Ја: 
| and for a vector function (СА) 
К як ме = fexto 
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48. We conclude our vector caleulus here. А much greater 
elaboration of details would certainly have been necessary, if it 
had been intended to present the subject with any degree of 

* completeness for purposes of practical aj 
certain obvious extensions moreover which suggest themselves 

й immediately from the work done bere in the foregoing pages j 

for example, 

(0) The discussion of the trilinear, and generally of the 
n-tinear symmetrical vector fuctions with a view to discovering 
the higher deri 

(ii) "The discussion of the improper intograls,—in particular 
integrals of functions having one or more infinities in the region 
of integration, and of finite functions through regions extending 
to infinity, and the differentials of such integrals ; and of Poisson's 
equation for vector functions. 

(ii) The consideration ої what Gibbs calls the determinant. 
of the linear vector function, and the development therefrom of 
the Jacobian and Hessian of vector functions. on the same lines 
that have been adopted here for developing the ideas of divergence 
| aud curl from what we have called the scalar and vector con- 
stants of the linear vector function. 2 














е operators. 




















we are always in the habit of thi 
(0 phy ideas that gave rise to them,—4o also form, quite apart 
. from their physical interpretations, the fundamental notions of 
the Abstract Calculus of Vectors, and supply us with a count 

` part of the differential co-efficient of a scalar function in the very 
‘real sense of giving us a basis of comparison of the rates of 
“change of the vector function from point to point of the field i 
` and further that this new mode of viewing them introduces con- 
‘siderable simplicity in the practical work of manipulation of 














| PART II pe 


‘Tue SrEADY Мотіох or a SOLID UNDER No FORCES 
IN LIQUID EXTENDING TO INFINITY. 






An attempt is mad this Part II to apply vector 
methods to the above problem. It ix just likely that it will 
be found to contain, especially towards the end, some new 
results which have not yet been worked out either with the 


cartesian or with vector caloulus, 








1. The origin О being fixed in the solid, if we denote by 
tho vectors R and G tho force and couple constituents of the 
“impulse” that would, at any instant, produce from rest the 
motion of the system consisting of the solid and the liquid, 
and by the vectors V and W the r and angular velocity 
components of the solid, then arguing as in Lamb's Hydro- 
dynamics, Chap, VI, $$ 119, 120, we have for the equations 
of the motion of system, 














AR 
RW ) a 
і and FRx VG xW =№ 


where Є, А represent the force and couple consi ituents of the à 
‘extraneous forces, the left hand sides being the rates of change 
of R, б when the ‘origin system’ (see Silberstein's Vectorial 
Mechanics, foot note, р. 69) bas the velocities V and W. 
If T is the kinetic energy of the system, we have 
2T=V.R+ W.G } = 
and further, R= v. T and G— V.T, 
` where у „Т and YT denote the gradients cf T regarded as a 0 
` function ої V alone and W alone respectively. [Lamb $ 122 


ea c ALT 
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3. We have to express first T in terms ої V and W and. 
then equations (2) will give us expressions for R and G in terms 


of V and W. 


If we put T, for the kinetic energy of the liquid motion 
alone and T, for the kinetic energy of the solid, T=T, +T, and 
we calculate T, and Т, separately 


3. To calculate T,—If U is the velocity potential ої the 
did. motion, we bave, if we take the density of the liquid. 
to be unity. 











2T, = fU U.da 
where the integration extends over the surface of the moving 
solid. 

Now the velocity potential U satisfies the following 
conditions :— 

(G) V*U=0 at all points of the liquid ; 

Gi) wv Uu. (V+W х) at any point P on the surface 


of tho solid, 7 denoting the vector OP, and я being a unit vector 
along the outward normal to the surface at P. —For the velocity 
of the liquid at the same point is VU, and the normal 
components of these two velocities must be the same. 

(iîi) yU=0 at infinity. 


OF these, condition (ii) shows that U must be linear in both 
Vand W. Itis also a scalar. Therefore it must be of the form 
F. V+F.'W, where F.F’ are two vector functions of the posi 
of a point (i.e, functions of r), but independent of V and W. 


Taking then U=F.V + ЕЛУ, condition (i) becomes 
уз (F.V)+ v ЧЕМ) =O. 
at all points of the liquid. Or, since V and W are constant 
vectors so far as the operation ої У? is concerned, 
VOy*FG WE, — [Se page 49.) 














ne 
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- Or, again, since V and W are perfectly arbitrary, we have - 
ViF-O and °F =O 2 (A) 


at all points of the liquid. 
Again F.V--F'W being written for U, our condition (ii) 
becomes м. 7(F.V + FW)». (V+W xr); 
or since У and W are arbitrary. — 
(В.У) =, y 


and mY (FW) иу xr 


at all points г on the surface of the solid. If for a moment _ 

we write $(êr) and ¥(8r) for àF and ЗК" respectively, we have 

by $30, р. 40, y(F.V)- &(V) and ЕЗУ) mud 

| therefore the above relations may be written NU 
- 





aN=nd(V)=Vid), | г 
and a. W x raa (W) e Wa(m)i 
or, a4 (m). в "a 
and гхи (н) н.у, } 
since У, W are arbitrary. 





Thirdly, condition (ii) becomes 
Vv(F.V)20 and g(F.W)-0 
ч or (У) =0 and y'(W)=0 E 
at infinity for all arbitrary vectors V and ХУ. ri 
© Conditions (A), (B) and (C) will 
vectors F, F' and so U being determined, we have 
? v re XT; fUv. da 
dt 5 i 
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4. To calculate T,, the kinetic energy of the solid ++ 
We have, ?T, = (V4 W x7). (V+W x 7) dm, dm being an 
`» element of mass of the solid at the point ғ and the integration 
extending through the mass of the moving solid, 
That із, 2T,=J[V.V+2V.W xr +W. (rx Wx7)] dm [see 
p. 6). 
=m V.V +@т у. Vx WWW), 
where denotes the mass of the solid, 7 the vector to it 
centre ої gravity (so that fr бта т), 
written for fr x (W xv)dm. Writing thi 
J(rrW —Wr.r)dm (p. б], we sce that « (W) is a self conjugate 
lineae vector function of W. Clearly also w(W) represents 
what the angular momentum of the solid about O would have 
been, if O had been fixed; and W.o(W)=constant, is, for 
variable W, the equation of the momental ellipsoid of the 
solid at O. 
5, We can now write down the expression for the kinetic 
energy of the system in terms of V and W. Thus, 
2T=2T,+2T, 


= JUV da + ДОМІ x r.da ну У 4-2mrN x W + W(W). 
For R and б, then, we calculate У.Т and Y,T from this 
expression for T. Noting that y(ar)=a and y(r.gr)=29r, 
if ¢ is self conjugate Гр. 28), and that V, U=F and V.U-—F', 
we write down almost immediately 

з Re2y.T-fUda- f FV.da fEW x гла +2mV + 2mW x 7 

and 








базу T= f F'V.da f Ur x дач FW x r.da + mF x V 
+2u(W). 

here are now certain relations among the integrals 

occurring in these expressions for R and G. Just to obtain 

these we | ч generally that if V is any constant vector, 





fev. (s y) pe^ N (da F), 


— — 
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Е, F' being any two vector functions, satisfying Y*F=0 and 
V*F'—0 at all points within the closed surface S over which tho 
integrations are performed. 

If we put, as before, 5 F=¢(8r) and 3F“ =y(ër), then da, E” 
=ф(4з) and Фа У F=¢(da). И then Cis any arbitrary constant. 
vector 


в 8 
CPV. (dag F’)=JF.CV.y(da) 
s 
=fF.Cy'(V)da 5 
в 
-[F.CV(F'V)da. 


T 

which by Green's Theorem=f7.[F-CY (F'.V)] dr, the volume 
integral being taken through the volume T onclosed by S. 
"Using now (i) $ 30, р. 38, and the relation vy*F'220, we have 


s T 
суку. (da. V) e еу (Е.С). v (FN) de; 


s T 
larly, C.fF'.V (Ja. уза [у СЕС) у CP Vir 
C being arbitrary, this proves our theorem. 
In the application of this theorem to our problem, we have 
to note that the region of integration would be that between 
the surface of the solid and a sphere of infinite radius, and the 
question of convergence of the integrals would arise, ‘This 
.—— question has been discussed by Leathem in “ Volume and Surface 
Integrals used in Physics," Scetions IX and XI. ‘ 

т. Using this theorem now and remembering our surface. 
condition (B), riz., фа) = da and y(da)=r x da, we bave, among 
` Elie integrals in the expressions for R and G in $ 5 У 





—e и © тетет 





اس 
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a ЈЕУ. da — fF'V.é(da) = [Р.Уф(дау= Р.У» x da, - 0 
Q0 [FW жена = JFW.rxda = ДЕУ. ¥ (da)= ЕЛУ (da) 





SSSR Wr x da -- 8) 
Hence those same expressions for R and G may be written 
R=fUda + mV + wW х= [Уда + mV]  [[E-Wda--mW x7], 
Gi [Ur x da + mF x У СМ = [f F. Vr x da + m7 x V] 
[IE Wr оуу). 
Wo write now R=p, V4 9$, W, where ф,, 4, аге the two 
linear vector functions, 
зу [E Na mV 
and Ф.М ПЕОМ ан xr. 


Since the conjugate of a. rû ix a г. û, (а) shews that фу ix 
> self conjugate. The conjugate of $, W 


=” W.da—wW x7 





ГРМ є x da 4 wr x W, by (у), 










so that ¢, V fF.Vrxda--wrx V. Hence we may write 
беф.У-о, У, 

where Ф.М = fF' Wr da +o(W). Since o(W) ін self conjugate 

(p. 78), (8) shows that $, is self conjugate. 
Thus we may write 

F R=¢,V+4,W } 
and G—4", V +4, W. 

where $u, $, are self conjugate functions. 


his fact alone,—of Фу» $, being self conjugate and Фф, and 
ў [И зетови oF conem be deduced directly from 


(3) 


Considering now the case where no extrancous forces are 
we have, putting £ A=0 in equations (1) of page 74. 
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ак 

аку, ж 

and 43 LRXVAGXW. v 


The three well known integrals follow immediately. 
бу Multiplying the first equation soalarly by R, we have ` 


R. Фе =(RRW)=0 
2 RR=constant. | 


Gi) 9. +в. - =(GRW) +(RRV)+(RGW)=0° 
5 







a R. 9 = constant. 


‘Phat is, the pitch of the wrench (R,G) which is EG is constant, 


Again, if re ERA, which we know is the perp. from O | 









on the axis of tho wrench (В.С), д 
d Rh ge x G+Rx AE] 


SS 





= gtx xo оона 





аа ет 
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Af then (47) denote the rate of change of r in spare, 






dr ж 
sp "ХУ +у=у— 


=V,, where V, denotes the com- 
ponent of V parallel to К, to the axis of. the wrench (R,G). 
It follows that th ія a fixed line in space, ite direetion 
being obviously constant from the first of the equations of 
motion. 





But there wax no special point deducing these results from 
the equations of motion, a» the fact they express, riz, the 
К сопмапеу of wrench (R,G) im case no extraneous forces aet, 
was obvious «f priori from the theorem that “ the ‘impulse’ of 
the motion (im Lord Kelvin's sense) at time £ differs from the 
impulse at time /, by the time-intogral of the extraneous forces 
` moting on the solid during the interval /—/,"—[Lamb $ 119)— 
` of which theorem it is only an an 1 expression that we have 
in the equations of motion. 























Gil) затенузом, 





2t aN AR y aw AG V 

М =. +7 У+О. р У 
у Ww AR ү tG 
T Using (3) р. 81, в ЛУ +6." reduces do V. A +W. G 





f. HR uw. dO e NRY 
a TE ev. +. лу = (VRW) AWRY) 


+(WGW)=0 

У UT constant, 
а ге integrals however are not sufficient to 
completely. We require three more scalar 


Б. А, 
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(or one vector) integrals for determining the two vectors V and 
W. The difficulty of finding them is avoided in two cases :— 


G) When R=0, in which case we can utilise the known 
solution of Poinsot in Rigid Dynamics. * 


(б) Steady motion, when A =0 and ~ =0, so that / 3 


ix got rid of altogether from our equations of motion. 


10. "The case R=0, is fully worked out in Lamb's Hydro- 
dynamics, $ 125. It may be interesting, just by the way, to 
put the solution in vector form, 





‘The equations of motion redace, when R =0, to 


яв у 
я таж. 


Again, putting R=0 in equations (3) p. 81, 
чо, М, 
G= ere WHW. 
Now ф, is self conjugate. Also, the conjugate of в ‘product" 
of linear vector functions being the ‘ product” of their conjugates 
taken in the opposite order (Gibbs Vector Analysis, Chap. У, 
p. 205), and further $, and therefore фу" also being elf 
conjugate, the conjugate of ф/,ф'Ф, is itself. "That is, G isa 
self conjugate linear vector function of W. і 
Hence, G4 V. (G.W) [жее р. 28.) pt 
Ix follows that Poinsot's solution is applicable. E 
11. But the саке of steady motion does not seem to hm 
received yet the attention it deserves. Two simple нан 


eases are well known—the permanent rst 
Г 





we have У 
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+ equations of motion of page 82 reduce to 


RxW=0 
ү Ü 
and Rx VG x W =0 i е 
2 ‘The first equation shows that R is parallel to W, or 
R=—<W, where sis a scalar. Substituting in the second, 


we have 
| GV +G)x W=0. 


Again, putting R= —zW in equations (3), p. 81, we have 
ap асо М, ХВ) 
"Were Wt М. 
W +26,W+0,W=0W, 
“Set Lt 0 md vete Ф - 
and 0 has been written for the linear vector function 1 
770, +0, FO 
We have, then, nWxWe0; e. (0) 





















and our conclusion is that in any ease of steady motion W should 

^de parallel to QW. 

19. By the rule for the conjugate of a product of linear , 
vector functions quoted in $ 10, wa soe immediately that 6,, 0,,0, 

` and therefore 0 also are self conjugate linear vector functions. 





— We recall now Hamilton's theorem of the latent cubie of a 
linear vector funetione—In general, for any linear vector 
functie б $, there are three vectors, say А,, Ass А,» of which the 
directions | are left unaltered by the operation of that function. 
Wage: 9a are the roots of the cubic equation РА 


= 5 
== 1, 
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— ww By у being any three arbitrary non-coplanar vectors, then, z 
ФА =A ФА, » and ФА, =9,А,. The cubic in g ix 
s roots the latent roots, and - 
the vectors A,,A,, A, which retain their old directions after the 
operation of 4 are called the axes ої ф. If ¢ ix self conjugate, й 
the latent roots are ull real, and the three axes real and mutually 
Perpendicular. [Killand and Tuit's Quaternions, Chap Х.]. 
Applying this theorem, we conclude that for any « there are 
three mutually perpendicular directions for W, correspon: 
any опе of which we may have a case of steady motion, 
Farther, if y is the Intent root of Q corresponding to the 
axis М, OW=yW, and we have iV+G=—yW, or 
| G-—v-yN. 


‘Thus our impulse i» given by 




















R=—-W 
e - 
and G2 — и i 
13. To construct a serew, therefore, such that the corres- ^ 


ponding motion of the solid may be steady, we find an axis 
of the linear vector func Ql) corresponding to any æ and 
“take W along this axis. V then is given by (5), р. 85 to 


ье олус. Draw the veotor WAV апа through. 


its extremity draw a line parallel to W. т line is the 

` axis of the screw. If then the motion is started by the 

impulse R= —-W and G= — V—yW, the solid will continue. 
to have the steady twist on our sorew, the angular velocity | 

















being W and the piteh Ду 






“It is easily seen the axis of any serew and the axis of the 
sponding wrench coincide. For these axis being pwralle 
and R respectively el. F 
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- drawn through the same point, жо that the lines coincide and 
either line is fixed in space. 








+ If y, р denotes the pitches of the screw and wreuch 
(R.G) respectively, 








14. We coneinde then that any = bei 
mutually perpendicular (though not necessa 
screws of steady motion are determinate, except only as to the 
magnitude of МУ ; and again that, by varying r, the directions 
of werews of all possible steady motions would be obtained by ` 
solving the vector equation П(г) xr=0, or (4*0, r4 10, r4- 0,7) 
xr=0, i 

^ funetion 0. 


specified, three 
intersecting) 














other words, by finding an axis of the linear vector 





It would be interesting now to enquire to what values of ж 
correspond the two cases mentioned in $ 11, ris, permanent. 
translation and steady motion with R=0, 


















Since, im any case, R= — МУ and W is not supposed to be 
infinite, R=0 would make «=0, О now reduces to 6,, 
G=—yW=-6,W, V= W 





Since, again, when the motion of the solid is one of 
translation only, МУ 20 and the screw. reduces to V only; and. 
since, in any case, R ix supposed to be finite, we see that in this 
case ж tends to infinity in such a manner that «W= 
‘The equations (8) of p. $6, being written 















Gow 6, Lo W) x R—0, 


© же nee that in this case 6, R x R=0, ie ф Rx R=0. 
Неше axis ої the wrench, and therefore the axis ої the 
screw also is parallel to r, where r satisfies dr x r=0. That is, 
parallel to an axis of 47 or to an axis of $., since Ф, and 

Pili ton az 
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These permanent translations, of course, could be worked 
out more directly by putting W=0 in (4) of р. 85 and (3) of 
p. 81. 


15. If we eliminate « from 





Vr буг 0,7) х0, 


we shall obtain the whole assemblage of lines, to one of which 
the axis of the screw corresponding to any case of steady motion 
must be parallel. 


Multiplying scalarly by 0,7 and 6,r respectively, we havo 





r.0, x r +0, rbr x r=0, 





and 278,78, xr +0, 


Та “r.a xr 
E 


or, 0,78 rx r) (вив x= 


rxr=0. 











rox", 


which is homogeneous and of the sixth degree in the tensor 
of r, and represents therefore а cone of the sixth order, to some 
generator of which the axis of every steady serew must be 
parallel. 


Obviously, the axes of the linear vector functions 6, and 
which we came across as giving the directions of screws for 
the special cases, $ 14,—all lie on this cone, for the equation їм 
identically satisfied if we put either 8,7 x rz0 or 6; x r0. 















be describing a helix about the fixed axis of the screw 
having the same pitch p (except of cotirse points on the 
which would move along the axis). For, if referred to 

origin O' on the axis, the position of any рой 


js specified by {OP =), the velocity of 


ses ділі 
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Hence, (W x г). (W x )e М xr. (Wx xt =Wxr, 


[W x(W x] =0 


h which shows that (W x 7).(W x) is constant, or the magnitude 
of W х г is constant. Thi » the distance PN of Р from the 
axis of W is constant. 











"The formula. for dr hows, moreover, that the velocity parallel 





to W is pW, and velocity in plane perpendicular to W is W xr, 
so that the motion in this plane is instantancously in a cirele 
(of which the centre is N and radius NP) with angular 








velocity i Hence P moves ín a helix of which the 
atan 
(00 piteh = =p. 
^ 
"A 17. For the maintenance of a motion of this type, appro- 









priate forces must be continuously acting on the solid, for we 
“know that the only motion a solid іх capable of under по forces 
is either one of uniform translation, or а uniform translation 
combined with а motion of rotation about a principal axis of 
the solid at the centre of gravity. We shall just verify that 
‘fluid pressures exert on the solid just the force and couple 
‘necessary for the maintenance of the sort of motion that we 
“ave bere. 





Considering the general case (where the motion is not 
cessarily steady), let £^, A” denote the foree and couple which 
uid pressures on the solid are equivalent fo. The linear 
alar moments of the solid are respectively 





(VW 7) and G, ат VW), зв) 
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these being just the terms of R and G that are obtained from Т, 
by the operation of V. and У. wepectively. For the solid 
alone, therefore, the equations of motion are 


| m 





=R, x W+? . 





and "> =R, xV +G WHA 


If we further pnt now R, =fUda and G, e fUr x da, жо have 
T R=R +R, G=G,+G,. From our equations of motion 
therefore of p. 32, we have 











яв, ^ 40 ЖОК Л. 
IR, eR x WE and “Os =R, x VG x WN 


Hence, fate +R xW and касі 
Rx VG XW. "s 


Considering steady motion now, for whic =0 and 


“Өз =0, we have 
2 (aR, x Wand MB, xV 40, x W, 
ог, using (4) р. 80. 
=k, x Wand =R; x V-0, x W. 
These formulae, for the special case of steady motion, could. 
“of courte be obtained directly by putting” =0 n =o 


Jin the equations of motion of the solid above. 
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1f, again, we refer the motion to the centre of gravity of - 
the body as origin, so that 7-20, we 








У бетухуУ 
and “=W x ((W), j 


where «КУУ) now i» the angular momentum of the solid about 
its centre of gravity. 


We notice that # and ХО both vanish only if either 
(î) W=0, or (i) У, МУ and «(W) are collinear. ‘The first ін 
the case of one of the three permanent translations. In the 
second ease, since W and e(W) have the same direction, W is 
along an axis of the linear vector function u, ie, along a 
principal axis of the solid at its centre of gravity. The axis 
of the serew therefore ix a priocipal axis of the solid at ite 
centro of gravity. ‘These two, of course, are precisely the cases 

` în which we expected 4 priori that б, ХО should van 

Wo show now that in the general case f, A’ of formula 

- (9) are just the force and couple that would make „the solid 

continue to have its werew motion represented by У, W at the 
‘contre of gravity. Since the motion of the centre of gravity and 

` the rotation of the solid about the axis W at its centre of gravity 

сап be considered independently, we show that the velocity V 

` of the C.G. is maintained by £', and then, regarding now the C.G. 
‘at rest, that the rotation W by itself is maintained by X 
4 comes to the same thing, that the mass-acccleration of the 
CG. ін equal to Є and that the rate of change of angular 
(momentum about еә. is equal to 3 

` For, if referred to a fixed origin O' on the fixed axis of the 
ми, the position of the C.G. is specified by ғ, then its velocity 
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домо =W x У, whieh is our Є. 






Again, the angular momentum of the solid about the С.б 
is (W), and in calculating its rate of change we take the C.G., 
which is now supposed to be at rest, as our origin, Since the 
“origin system” [see Silberstein's foot note cited on p. 74) ін 
fixed in the body and rotates with it with angular velocity W, 





the rate of change in question = Ets (м) =W 


| =Wxu(W), (> у =0), : 
which is our A". ‘ 
1 

le 18. Having thus considered the general character of steady 










we would next turn our 










to summarise here a few propositions on the linear veotor funetioi 
which we shall presently bave occasion to use. ‘The more ` 
important ones are taken directly from Kellaud and "Tait's 
Quaternions, Chap. X. P 
(i) From the well-known relation - 
козу) = (By + (тед + (сайуу, 
it follows easily enough, — say, by writing 4(7) in the form 
Ar + a.rb +e xr of page 24—that 
(B) = (BY) Ha) + Cr) (8) + Cr )6)- 
“Hero a,f,y are any three arbitrary non-coplanar vectors, and 
К апу fourth vector, and ф denotes а linear vector function. 
ER explained on p. 6 staud for а xy, ete- 
Gi) For inverting the function s) we have, if wed n 
Function by ¢ адр 
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where (ay) = ($2986), and A, м апу two vectors. 
If we introduce the function $ by the definition, 

Оха) Ao HH, 





we may write symbolically пф" =v. 





(їн) For inverting фу + gr, ог ($ + 7)r, where g is a constant 
scalar, we have 


Отто m"g* 07) (+d О н) = (9 9) (9 + дун 
E =(+9x+9") Axr), 
where m, ют", m” have the values defined on p. 86, wir. 


Ву) = (ф°ФЙФу), т (а/8у) = х(9ф8%у), m" (afr) = Х(Вуфа), and 
` the function x is defined by 


ХА х p) m 4 x и АХ фо 





Me 






(iv) Integrating dr over the surface of а parallelopiped of 
С whieh the edges are the vectors o,f we get easily enough 
С perdam (Вун). 


Hence our "= =realar constant of $ 


(v) We have also for any vector A. 
m" hd + АА, 
symbolically, wx. 

(Wi) We already enunciated Hamilton's theorem of the 
Intent cubic on page 86. We only note here that if for any 
function ф, m=0, the product of the three roots 7,» 9, and у, 

‘the latent cubic vanishes and therefore one of the roots, say 
š o. Tt follows that for such a function there exists a 
, such that 43, =0. Conversely, if for a function ¢ we 
vector A, such that ФА, =0, т for that function must 
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Keeping to Gibbs’ notation, we shall speak of m for any 
function $ as its determinant and denote it by | Ф | 

(vii) About the determinants of linear vector functions, 
we have the theorems that the determinant of a ‘product’ of 
linear vector functions is the product of their determinants, and 
that the determinant of tho ‘quotient’ of two linear vector 
functions is the quotient of their determinants. Thus the 


dah le 


ме the determinants of ¢,, aU ss ee [See Gibbs, 

р. 2.] 
| (vii) The inverse of the ‘product’ of any number of 
linear vector functions is the ‘ product’ of their inverses taken 
in the opposite order. (Gibbs, p.293]. Thus the inverse 
Of dideds=(dideds) 6; 7 24 

10. To examine now the stability of any particular mode. 

of steady motion, we consider ая usual the effect of n small | 
disturbance given to the system, У, W being the linear | 
angular velocities of the solid for the steady motion, and В, 
the corresponding impulse, if ЗУ, SW are the additional linear 
| and angular velocities imparted to the wolid by the. disturban 

the total impulse that would generate from rest-the new 
disturbed motion is represented by Ae 


RR, (V КУ) ФУ +8W) 
and G+3G=¢',(V +8V) + Ф,(%У +8W), 

no that зва (FV) + 4, (8W) 
and 6G=¢',(5V) +A, (W). 3 
Our equations ої motion of p. 8? now are 


f determinant of ф,ф,ф, 





sif mmm. 





































Фе. 


© and (6 +26) =(8 +88) онуно 
Resim 
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which, by the condition of steady motion ($ 11, р. 85], 
reduce to 


M ФК. 
GE Rx êW вх, 


and ЖӨ (КВУ +R x V НО x 3W 456 x W, 


if we neglect the terms 8Rx5W, 3R х 5У and 56 х 8W, 
Using (7) of p. 87, these equations farther reduce to 





GR SOR +8) хуу, я 

(10) 

and O — (GR 18W) x У + (8G + 18V + уйуу. 

(30. Let us put now ЗУ ет" and ВУ еко, where v, 10 

wo vectors independent of /, so that 
3R=R'’e"' and 6G=G'er', 

where Ri=d,0+4,0 


ho р b „взу 
and Офф: 














ма HB uer Ото 
E (Equations (10) then are identically satisfied provided 





, һВ'=(В/'+ с) x W, an 
and со В РОБИ zs 
autem 
- (14) 
і Z ; 
беу tex M 


To n EIC 
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From (13) now we bave D азу 


and —=Ф,®р'Ф,е+Ф,е=0 . ae 






It follows that the latent cubic of the fanction—#,@;'@, ФФ, 
has a zero root and w is the corresponding axis. Thus the 
determinant of the function vanishes and we have 








о) [8,8—9,0,6,8) x [6,5—,9; 9, y] =0, 
[See (її) and (vi), $ 18] 
а, В, у being any set of three non-coplanar vectors ; that 
хе, у 
у}+Җ[Ф%,«Ф,Ф;'Ф,8 хФ,Ф'Ф,у]=0, 
Dividing out by (a8), and noting that by 5 18, (vii) p. 99, the 















determinant of eem e ا‎ Sore Eure: 





f 
| 18, | - LL Pel ру < tst] 
\ hh ха „Вх®, yJ=0 .. (17у 


This relation (17) gives us an equation for the appropriate 
values of н, and only when the roots are all imaginary, would 
the motion be stable. 

91. Summary of results to be ргогей--Ль iw very tedions 
mow to calculate all the co-efficients of the several powers of т 
їв thin equation, and even when these co-eiicionts aro calculated, | 
it is found impossible to come to any definite conclusion as to 
‘the conditions of stability. The solution of the problem with 





З 
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for is of the seventh degree, and that the co-efficient of я" ix 

determined solely by the nature of the solid and does not vanish, 
+, unless some special limitation is imposed on its shape, вау, by 

way of symmetry. If will follow that, being of an odd 
degree, this equation must have at least one real root, and that 
in general therefore, for any arbitrary г, we shall not have any 
stable steady motion at all, We shall next show that for some 
special values of +, the term independent of » in our equation 
vanishos, and then neglecting the root #=0,— of which the 
offect is simply to add a constant vector each to the values of 
8V and 8W, and which therefore does not affect the question of 
stability either way,—we shall have our equation reduced to 
one of the sixth degree; and since an equation of an even 
degree may bave all its roots imaginary, it is possible now for 
the motion to be stable. Instead of groping about then for the 
stable sorews among the infi system of steady sorews, we 
shall be sure of this one fact that if any stable screws ate there, 
they would be found only among the group determined by 
those values of г which make the term independent of я 
vanish, It will be shown again that in general there are six 
such values of s, and our conclusion would be that in any саке 
there cannot be more than 18 (—6x5) stable steady sorews. 
We shall show finally that «= is one of there віх values of г, 
and we shall conclude by working out this case completely. 

33. To obtain the degree in » of the loft band side of (17), 
let us first express R’ iu terms of # from the first. of the equa- 
tions (12) of page 101. 

Writing фВ/ for «R'—R' x W for the time being, this equation 
is éR'znR'—R'xW-aexW. Hence by (i) of $18, р. 97, 
we haye 








































mR'— атф (wx У) афо x ФУУ, 
Where фло ix the conjugate of due and is therefore nur Fie x W 
ы #WanW, and 


j) = (чаа ж W).(ng—BX W) x (нуу x W), which, 
7 itten out in pu [1 жщ Apa мые 
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‘Thus we have, n(n* +W.W)R'=r(me+wx W) x nW, ї 
or, (at + W-W)R'—s[n(we x W) +(e x W) x W] " 
вне x W)4-WaeW — W.Ww], p. 6, ^ 08) 5 5 


23. Now from (11) p. 101, R'=¢,0+¢,, Therefore, 
TP ФЗ". Kquating this value of v and the value 
in (15), р. 101, we have, Бу using (18), | 





wee; Гафт (u ж W) 4 Wed W 





we hm TEW 

улут) 

Te (нефт de ну (e x W) 
+ WWE ba ee We tW] 

Looking wp the form of 4, now in (14), p 101, we ме that we 

can write зе WW (n^ Ant РР РР where 


A, B, C, D denote certain linear vector functions which do mot 
involve n in their constitution. We have, for instance. 


еф ds bees * a 09) 

















nd again, it is easily seen also, that 
—Юю=ф,[%# МУФТ (9, не У ву  W]x V 
E ну Wes (b, + У ФУ 
=үлї{(Ф.+)ехУ e (s FOOT Gs + е W] 
—ум{—У+(Ф',+)фФ: W1x W 
зум. + е x V + ($, E) W) 
ure xW, A 3 
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24. Our X[@,0.0,8%,07'4,7] of (17), р. 102, then, is of 
the form 


a т ж [WAH B')a (84^ + B^) x (nA n* B+ -D)y]. 
M where nA'4-B' has been written for ®,, чо that 
q А'а=ф,а, and ~Ba=(p, +«)ax V ($, куда x W 
9 Obviously A’, B' also do not involven. We may write therefore 





Ge By) Met. Bx ter py] 














рум [e aunt ba n! Erat ea neu] 


whore tho a's are all independent of я. The values of a, and ay 
may be written down immediately. Thus 

(oBy)a, —3[A'-A'B x Ay] — 3|, вх M de ey c (21) 
and (ауа, зад В'є BB x Dy] ^ (0) 


35. Again, with the same notation, 





vB x. 





sy) 





сю 3 (2,5. Ф, 


BI en S[(mA'4+ Ba. (n* A--n*B--nC-- D) 
: x(t A+  BA-nC -D)y], 


=3. ABx Ау оліфа ө. 
xe. t ey oe (28) 
(9By)a', —X[B's. DB x Dy] -- (24) 


pe 
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26. We try to form an idea now of the other terma i 
\Ф, |, 19, | and |, | that occur in (17), р, 102, | 


@) We note in tho first instance that ¢(r) being any liner. 
vector function and a any constant vector, the determinant of 
i фт xa vanishes, For, 


і (фаха). ($B x a) x (фу x a) (фах а). а(фДфув)=0, [page 6.) 


It follows that the parts independent of n in ФТу and Фут, wis, 
фуу and—($,+2)rxW have their determinants equal to 
zero. 


vw 








| (йу We note again that if $,, $, are any two linear vector 
[ functions, the determinant of nġ,r+¢,r is easily found by writing 
| out the expanded form of 













суу ОНА), ФВ) х (nd. MIs 


to be m,n? +n" tu'n tmy, where m,m, nro. the determinants of 
" Фу, $, respectively and p, м' are determined from 


n(oBy) з 3Гфуа. d$. Bx Фу) 
and нау) =21%,9. $, 8x $,yl 


It follows that the terms independent of n vanish in both 
|4, | and |, |. We find, in fact, by working out 









Quy nte he WD. Fp, EX) n nam hax WO 


as in § 22, that 
Те, | em, пин ЗУУ), 
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where m, is the determinant of either of the two conjugate 






















functions ф,, 4^, and m, is the determinant of фу; and the Ws, c's 
+ and d's are all independent of mand in general none of them 
vanishes, 
27. For our equation in т then, corresponding to (17) р, 102, 
we have 
i man" bunt cun d, 
( (mam? ban? e nd, Стун" bun eu) 
E m WW ° 
2 mant. ta, 
| www - 
М ni +... +a", 
| ж ер МУ 
=0. ^ x = - (85) 


‘Multiplying out by erm we get an equation of tho 
seventh degree in n. 


"The co-efficient of м 
mato, Mi 


=1%1- ERAI 


зару Зі ati eX! ba 

і -Is 

_ by (21) p. 106 and (23) p. 107, and (vii), p. 99. 

1 peer |, $ 26, we seo that this co-effi 

3 $, | =— | ô, |ô, having same definition 
o 795 Obviously, our co-efficient of n” does not vanish, 

th пк $i dy an $, are given in special ways, 

that is to say, some special restriction is imposed on 
of the solid, for these functions are solely 

d by the form of the bounding surface of the solid, 
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28. We may write down also the (erm independent -of a in 
the same equation. ‘This term directly 


=(W.W)"d,—2, W.W a^, — WW EX D 
Now if, as in 5 24, we write n АВ for Ф, d, by (ii) § 28 in 
he determinant of B. That is, 
4 (o) - CB'aB'BB'y). 
The values of a, a’, are given by (22) and (24), vis., 
Саву) o. — S[ Da. BB x Dy] =W.W3[B'a BB x D'y). 
(ufly)o" ,  X( В'». DB х Dy] (WW)? XL Da, D'A x Dy), 
if a new function 1У ix introduced by the definition 





—D'r= үзүр Ре which by (20) p. 108 
mb +), к V +[(ф, HY WF] x We 


(240, W +0,У) 





where | Рос 
Henco, (W.W)'d, —a, W.W ка", 


=" oremus saremo y) авер) 





=(W.W)* | B'—D' | + | Ю"1 J, by Gi), § 26. 
‘But (B.—D/)r- [(0—y)r4- WF] x W, 
E ues eie 1218—17 | =0. 





у STEADY MOTION OF А SOLID 85 





| 
| 39. For further simplification of this term, we calculate 
[| gonerally the determinant of 


| E Pax Vy W, 
b Фу, ф, being any two linear vector functions. 


Since the determinants of both $,yx V and $, x W are sero 
\ by (1) $ 26, the doterminant in question ix by (ii) 5 26, 


= a (фах V). вх V)x (huy W) 
жах У), (4, 8x W) x Guyx МО) 


(0 Now X(4 ax V) x ($48 x УМ (doy W) 
фу оф, AVIV. (фу х W) 













Е 
ту =W ху. HOT Visas Bey). by б 
c =з ху. MeV mi Cay), m, 
С Agnin, X($, «x V). (4,A x W) x (be W) 
1 зе V) Ми) А 
=V XW. Зв, Ве) Я 
уху. өүө Wib BD 

=V xW. ф,ф;'ЎЎ т, (98у), m, being `$, | 





Tho determinant of $,rx V+ фт W is therefore 
2l = ХУ. [m dedit Vom $i W] s dH) 
na SW V. (4,4, Vf, WW], in the notation of. ens 


thie result now, we easily vnleulate d, and 
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Thus, if we write nA" +B" for ®,, 4, is by (il), 5 26 the deter- 
minant of B”. That is, 
—d,=determinant of $,r X V ($^, аж W 


=WxV. [HL tO, VOY Fee, + )W], by 
(iii) § 18, p. 97 


Again,— | D' | determinant of (фу + z)r x V+frx W, 

(020), p. 110) 

if fris written for the linear vector function ($, +0)r+W.rF. 
Hence, wo have— | D' | =W x V. [f(Y, +4 x, +4*)V 


(Фи), 
where g denotes the y function for /, that is, po 


захв) = Лау = (бака) a+ Wi Fa] x ($40) в 
+W. FA]. 


31. For the expression (27), с, only remains to be calou- 
lated. This i» done directly. Thus by (ії), $. 26, : 


с, (обу) = 2ф,о. [(Ф, + =) Bx W)) x [Gh + =) yx МУ) 
= Fhe. W (0+) BX (фу + y. W] 1 
= Sua W [Gy HW.) By]: W, [see С, p 97] _ 






which by (i), $18 = (98у) W. [а%%, W нива) 
(^ since by (ii), $18, dYa W = ms We 
2 emt м. = WoW. fans W 
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where сз, d, and | D'| have the values calculated in the last 
two articles. As explained in $21, it is necessary for stability 
that this term should vanish. A necessary, — uf course, not 
sufficient — condition therefore that any one of the threo steady 
(^. sorows corresponding to = should be stable, is that є should 


satisfy w.wqp 0 =0 КУ м (28) 
a 





Since now z occurs in the first degree in У [(5), р. 85], in 
the second degree in Q Гр. 86] and F [(26), p. 110] and in the 
fourth degree in 7, it appears from an inspection of the values 
f of cy, ds and | D’| given in the last two articles, that (28) 
> represents an equation of the sixth degree in s, and will, iw 
general, therefore determine only six values of г. No » other 
than these six can possibly determine a stable steady serow. 


33. We prove now that «=0 is always one of this set of 
six values of х. It їв only necessary to show that (28) is satis- 
fied when «=0. 


For, when 420, the expression for —D'r in (26), p. 110 
reduces to 















dur XV FP ST Sx W, 


Е being zero, and ¢ +0 reducing to $s + 
in this case. Honce now by (iid, p. 112, 


| = |D | = У. Ст, 4.47 $67 V— | во, | 
з HOTS) Ww). 





0. to Pab Pa 





| But, by (vii) and (eii 





$15, 


Фе 





у живо» | = TE and Суто 
+ 
» evo zn зи 


= Me Wx V. (hs V e by (i), $18. 


—Ó 
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Also, putting == 0 in the values of ep, 4, in $630, 31, we 
have now с, =m, W. W 


and —4,—W x V. (97, Уф М) 


It follows that W. W | D' | — 26 





* —0 identically when 





2-0. 


34. Thus the case »=0 satisfies our necessary condition of 
stability. It is a relevant enquiry thon if any one of three 
serews corresponding to this сало is stable. This is the case 
($14) when the impulse producing the original ру motion 
reduces to a couple alono, G = —0, W =—yW ; R=0. 








For equations (12) of p. 101, we write now 
ноз вих М 
«бе вх V+ (G'+yo) xW. 
From the first, В' = 0 (for no vector сап be perpendicular to 
itself), and therefore the second reduces to 
nG'=(G' + уе) х W am .. (09) 
Putting again R' 
vm фуфло, О = — dite hu = 





in (11) р. 101, we have 








Writing 4 for 6, for convenience, we have G'= 
and n now (Фу ex W =0 
er nite (Фо y)m x W =0. e (30) 


Putting the determinant of the function on the left band 
sido of (30) to zero, we shall have oar equation for the appro- 
priate values of з in this case. 

The determinant in question ie, by (i) $ 26 
т т? bn? +en +d, 
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where, —/= determinant of ( —y) «x W =0, by (г) $ 265 
== determinant of $ = | ,„ | 
I W. GW —yW. ФУУУ +mW. W, just as in the 
calculation of є, in $ 31; 
and —4 (ау) = (фа x $8). (9—0) yx W], which we find to 
be zero when we write it out and remember also that 4 is self 
conjugate ; so that 4—0. 
Hence, our equation for » now is 
maden W. ФУУ УУУ. ФУУ + Уу, W]20. .. GI) 


[We may note, by the way, here that the co-efficient of м? 
in this equation, m or | „|, is the same (but for the 
sign) as the co-efficient of wî in the equation for м in the general 
oase (нео $ 27, p. 109). ce this last co-efficient in any case 
is independent of г, and since equation (31) ix only the reduced 
form, when 2 =0, of the general seventh degree equation, ме 
see that when =0, not only does the term independent of м 
vanish in the general equation, (as we bave proved in $ 33), 
but the co-efficient» of all terms from #* downwards vani 
too}. 

35. Now one root of equation (31) ік #=0, and the other 
























two are wiven by 
wt — 1 (PW. ФУУ У + mW, W) 
since ФУУ зо МУ, and by (V) $ 15, xz m" — o, this simplifies to 


wt = — UM [age терн e). 





‘The motion їх stable, therefore, if y i» so chosen that 
ر‎ [249—w"y*-m] . e. (38) 

je positive, Now, by Hamilton's theorem ої the latent cubie, 
($ 12, p. 86) у may have any one of three values, ríz, the roots, 


say уз, Уз, Yar OF the equation 
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ТИ БУУ. AI ч my rur 
, у equal t any;one of Phere Fols e» s. 
(32) besan ; Ў 


riori Gn coser) nad 


ЛЕ утар 


ти) Gur) ? 7 (взу 


fln ately „егуз, therctore, parallel to that axis of, ої 
to Ж latent root y ie table if, thie expresion (33) 


a R being zero, the enenzy ої the steady motion for 
wy d GW = — py W. W, which must in any само be 
positive, ft follows that all Diem 7 The oo 

| therefore) їх positive, itty, ix planes "the greatest 


б ive, 


